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SUBMANIFOLDS THAT ARE LEVEL SETS OF SOLUTIONS TO 
A SECOND-ORDER ELLIPTIC PDE 

ALBERTO ENCISO AND DANIEL PERALTA-SALAS 



Abstract. Motivated by a question of Rubel, we consider the problem of 

characterizing which noncompact hypersurfaces in R" can be regular level 

sets of a harmonic function modulo a C°° diffeomorphism, as well as certain 

generalizations to other PDEs. We prove a versatile sufficient condition that 

^^ ' shows, in particular, that any (possibly disconnected) algebraic noncompact 

^SJ , hypersurface can be transformed onto a union of components of the zero set of a 

harmonic function via a diffeomorphism of R". The technique we use, which is 

a significant improvement of the basic strategy we recently applied to construct 

h^ ' solutions to the Euler equation with knotted stream lines (Ann. of Math., in 

■^^ , press), combines robust but not explicit local constructions with appropriate 

• 1 global approximation theorems. In view of applications to a problem of Berry 

i-i^ ' and Dennis, intersections of level sets are also studied. 

-)— > 



^ . 1. Introduction 

> 

A long-standing open problem on level sets of solutions to elliptic PDEs, formu- 

00 . lated by L.A. Rubel and included in [H Problem 3.20] and [ITl Problem R.7], is 

to characterize harmonic functions in K" up to homeomorphism. That is, given a 
continuous function / on K", one would like to know whether there is a homeomor- 



f^ I phism $ of R" such that u{x) :— /($(a;)) is harmonic. In particular, a necessary 

^^ ■ condition for the existence of $ is that the collection of all level sets of / must be 



homeomorphic to that of a harmonic function. 

Actually, Rubel's problem was motivated by the foundational results of Ka- 
plan [16] and Boothby [6] in the theory of foliations, who studied a less stringent 
l^% ' version of Rubel's problem in the case of dimension n = 2. This relaxed formulation, 

S . which is more natural from a topological viewpoint and goes back to Morse 1191, 

is tantamount to characterizing the collection of level sets of a harmonic function 
up to homeomorphism. Equivalently, given a continuous function / on R", the 
question is whether there exist a harmonic function u in M" and a homeomorphism 
$ of R" which maps connected components of any level set /^^(c) into connected 
components of level sets of u. 

Even with Kaplan and Boothby's relaxed formulation, Rubel's problem is wide 
open. The easiest case is that of dimension n = 2, where complex analytic meth- 
ods are of great help in the study of the singular foliations defined by harmonic 
functions [121 M HZ]- Particularly, these authors used Stoilow's theorem to show 
that, for any continuous function / on R^ whose level sets satisfy certain neces- 
sary conditions, there are a homeomorphism $ mapping the unit disk D^ onto R^ 
and a harmonic function v on D^ such that v{x) = f{^{x)). In regard to har- 
monic functions on the whole plane instead of on the disk, these authors derived 
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Figure 1. From the left to the right, an infinite jungle gym, a 
torus of infinite genus and a torus of genus 5 = 3 and N = 3 ends. 



nontrivial topological obstructions (both for the relaxed and Rubel's formulations 
of the problem) on the admissible families of level sets of / as a consequence of 
Picard's theorem, and constructed infinitely many non-homeomorphic families of 
level sets of harmonic functions. A sufficient criterion for / to be homeomorphic to 
a harmonic polynomial was obtained by Shiota |24| . 

To the best of our knowledge, there are no results related to Rubel's problem 
in any dimension n ^ 3. The reason for this is twofold. Analytically, all existing 
results in the case n = 2 have been obtained using complex analytic techniques, 
which are no longer applicable when n ^ 3. Topologically, the situation is much 
simpler when n — 2 because any connected component of a regular level set of a 
function u on M^ is homeomorphic to a line or to a circle, and the latter possibility 
must be obviously excluded in the case where u is harmonic. On the contrary, in 
higher dimension there are infinitely many topological types of level sets and it is 
not obvious which ones can actually be level sets of a harmonic function. 

Our objective in this paper is to study which hypersurfaces L in M" can be 
regular level sets of a harmonic function up to homeomorphism. In particular, we 
aim to explore to which extent there is a large collection of hypersurfaces that can 
be level sets of a harmonic function. Roughly speaking, this would measure the 
complexity of the (singular) foliations one would need to consider in any approach 
to the higher dimensional Rubel problem. Other than the obvious obstruction that 
the hypersurface L must be noncompact and the analysis of explicit examples of 
harmonic functions, there are no results on the admissible topological types of L. 
For instance, in dimension n = 3 it is easy to construct explicit harmonic functions 
having level sets with a connected component homeomorphic to elementary surfaces 
such as the plane or the cylinder, but nothing is known regarding more complicated 
objects as in the following 

Question 1. Is there a harmonic function u in M^ such that m"^(0) has a connected 
component homeomorphic to the genus g torus with N ends, to the torus of infinite 
genus or to the infinite jungle gym? (cf. Figure]^. 

Our main theorem provides a sufficient condition for a hypersurface to be a 
connected component of the zero set of a harmonic function in R" up to a diffeo- 
morphism, thereby furnishing a positive answer to Question [T] as a particular case. 
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More generally, we will prove this result for the equation (A — k'^)u ~ 0, with k 
a real constant, and for possibly disconnected hypersurfaces. A somewhat related 
result, which gives sufficient conditions for a set of curves in the plane to be home- 
omorphic to the zero set of a harmonic polynomial, has been recently established 
in [10]. 

Thinking for concreteness in the three-dimensional situation, the strategy of the 
proof is to construct a local harmonic function having a level set diffeomorphic to 
the surface (for instance, the genus g torus with N ends, the torus of infinite genus 
or the infinite jungle gym, as in Question [IJ whose gradient and nearby level sets 
satisfy certain (rather subtle) geometric conditions. Once this function has been 
constructed, we show that there is a global harmonic function that approximates the 
local solution uniformly in a neighborhood of the surface, and prove that this global 
solution has a level set diffeomorphic to the surface using a noncompact stability 
theorem and the geometric properties of the local solution. In order to carry out 
this strategy, it is crucial to deform the surface to ensure that its geometry is 
controlled at infinity, for example by transforming the collared ends of the surface 
into straight cylinders (as in the case of the torus with N ends) or by making 
the surface invariant under a discrete group of isometries (as in the case of the 
infinite jungle gym or the torus of infinite genus). As a consequence of this, the 
homeomorphism that appears in Question [T] will be obtained as the composition of 
a 'large' diffeomorphism used to control the geometry at infinity of the surface and a 
'small' diffeomorphism related to the approximation of the local harmonic function. 
A simpler application of the philosophy of combining a robust local construction 
with a suitable global approximation theorem has been recently applied to study 
periodic trajectories of solutions to the Euler equation [^ . 

Because of this two-stage construction of the diffeomorphism, the theorem can 
be most conveniently formulated in terms of a natural n-dimensional analog of these 
surfaces with controlled geometry at infinity, which we call tentacled hypersurfaces. 
The precise definition of the latter is given in Sections [5] and [51 For the purposes 
of this Introduction the reader can stick to the three-dimensional case as above, 
keeping in mind that, for instance, any hypersurface obtained from a compact 
manifold by removing patches and attaching 'nicely embedded' collared ends (the 
'tentacles') is diffeomorphic to a tentacled hypersurface. 

Theorem 1.1. Let L be a tentacled hypersurface o/R" (possibly disconnected and of 
infinite type) and let k be a real constant. Then one can transform the hypersurface 
L by a smooth diffeomorphism <I> o/M", arbitrarily close to the identity in the C^ 
norm, so that $(£) is a union of connected components of a level set u^^{cq), where 
the function u satisfies the equation (A — fc^)u = in M". 

It should be emphasized that the class of hypersurfaces that can be realized as 
a tentacled hypersurface modulo diffeomorphism is deceivingly wide; in particular, 
this includes all real algebraic noncompact hypersurfaces, which allows us to state 
the following 

Corollary 1.2. Let k be a real constant. Any (possibly disconnected) noncompact 
real algebraic hypersurface L of K" can be transformed by a smooth diffeomorphism 
$ o/M", not necessarily small, so that '^{L) is a union of connected components of 
a level set of a function u satisfying the equation (A — k'^)u — Q in M". 
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The subtleties that appear in the proof of these resuhs are related with the 
necessity of combining the rigid, quantitative methods of analysis with the flexible, 
qualitative techniques of differential topology. A simple yet illustrative example is 
the necessity of allowing for a small diffeomorphism $ 7^ id in the above theorems: 
indeed, it is well known that the curves X2 = x'l, which are all diffeomorphic for 
any real s ^ 1, can be a connected component of a level set of a harmonic function 
in M2 if and only if s = 1 or s = 2 [TT] . 

After discussing the level sets of a harmonic function, we will next consider 
regular joint level sets of m harmonic functions, i.e., transverse intersections of their 
individual level sets. Unlike a single level set, these intersections can be compact 
provided the number of functions m is at least 2. Again, nothing is known on 
this problem, particularly if one is interested in sets diffeomorphic to pathological 
objects as in the following 

Question 2. Let L be an exotic sphere of dimension m, so that L is a smooth 
m,anifold homeomorphic to the standard m-sphere but not diffeomorphic to it. It 
is well known that L can be smoothly embedded in R^™. Are there m harmonic 
functions Uj. in R^'" such that the exotic sphere L is diffeomorphic to a component 
of the joint level set wf ^(0) n • ■ • n u;;^^{0) ? 

Our main motivation to consider intersections of level sets comes from Berry 
and Dennis' question [5] of whether any finite link (that is, a collection of pairwise 
disjoint knots) can be a a union of connected components of the zero set of a 
complex- valued solution to the PDE {A + k^)u — in R^. Their question, which is 
of interest in the study of dislocation structures, is a particular case of the following 

Question 3. Let L be a locally finite link in M.^ . Are there two solutions Ur of the 
equation (A + k'^)ur — Q inM? such that L is diffeomorphic to a union of connected 
components of the joint level set u^ (0) n U2 (0) ? 

We shall next state our main results on joint level sets of solutions to an elliptic 
PDE, which in particular provide affirmative answers to Questions [2] and [3l Here, 
our results and the method of proof depend greatly on whether the joint level set we 
want to prescribe is compact (or, more generally, a union of compact components) 
or not. The case of compact joint level sets is less involved and can be tackled for 
a very wide class of elliptic equations: 

Theorem 1.3. Given an integer m^ 2, let L be a locally finite union of pairwise 
disjoint, compact, codimension-m submanifolds of M" with trivial normal bundle. 
Take m real constants Cr and m second-order elliptic differential operators Tr with 
real analytic coefficients. Then one can transform the submanifold L by a smooth 
diffeomorphism $ o/M", arbitrarily close to the identity in any C^ norm, so that 
<i>(L) is a union of connected components of the joint level setui (ci)n- • -nu^^^lcm) 
of functions Ur that satisfy the equations T^Ur = in R". 

Our approach to the case where the prescribed joint level set has a noncompact 
component strongly depends on the method of proof of Theorem II. 1[ so the result 
will be stated in terms of the codimension-m analog of a tentacled hypersurface, as 
defined in Section [T] To this noncompact submanifold we can also add a union of 
compact components of the kind dealt with in the previous theorem: 
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Theorem 1.4. Given an integer m ^ 2, let L he a tentacled submanifold of codi- 
mension m in M" (possibly disconnected and of infinite type) and let L be a finite 
union of pairwise disjoint, compact, codimension m submanifolds ofW^ with trivial 
normal bundle. We assume that L and L are disjoint and fix m real constants k^. 
Then the codimension-m submanifold LU L can be transformed by a smooth diffeo- 
morphism $ o/R", arbitrarily close to the identity in the C^ norm, so that ^(LUL) 
is a union of connected components of the joint level set MJ~ (ci) ■ ■ ■ D u~^^{cm), 
each function Ur satisfying the equation (A — k'^)ur = m M". 

The proof of these results also rely on the combination of a robust local con- 
struction with a suitable global approximation. It is worth emphasizing that our 
treatment of noncompact joint level sets is considerably simplified by the fact that 
the diffeomorphisms that appear in Theorem II . 1 1 can be chosen C^-small (and not 
only C"-small). In proving this bound for the diffeomorphism, we make essential 
use of fine C^-estimates for some Green's functions; had we restricted ourselves to 
the standard (C^) gradient estimates, the diffeomorphisms would have been only 
C''-small and the treatment of joint level sets would have been much more in- 
volved because the transversality of the intersection of hypersurfaces is generally 
not preserved when one applies a C^-small diffeomorphism to each hypersurface. 
The lack of appropriate estimates for higher derivatives of the Green's functions 
also lies at the heart of why, in the noncompact case, we can make all the above 
diffeomorphisms C^-small but not C-small for all p. 

The paper is organized as follows. To begin with, in Section [5] we provide a 
guide to the demonstration, where we present the global picture in a less technical 
manner and explain the relationships between the different parts of the proof. In 
Section [3] we prove the novel noncompact stability theorems that extend Thorn's 
isotopy theorem to the noncompact case. These results play a key role in the proof 
of Theorems 11.11 and 11.41 snd are proved using a semi-explicit construction and an 
appropriate gradient condition to deal with the lack of compactness. In Section 2] we 
prove the necessary global approximation results that are needed in the rest of the 
paper, which generalize theorems of Gauthier, Goldstein and Ow [12] and Bagby [3]. 
Once these crucial preliminary results have been established, we present the proofs 
of Theorem 11.11 and Corollary 11.21 in Sections [5] and \6\ and those of Theorems 11.31 
and 11.41 in Section [T] For completeness, in the Appendix |X] we consider the case 
of a single compact level set: while this situation cannot happen e.g. in the case 
of harmonic functions, which is the main thrust of the article, we show that any 
compact hypersurface is diffeomorphic to a level set of a solution of the equation 
(A — q)u = in R" under fairly general hypotheses on the function q. 



2. Guide to the paper and strategy of proof 

In this section we shall give the global picture of the proof of the main results, 
which will enable us to show the connections between the different parts of the 
proof without technicalities. We will also present the organization of the article, 
explaining the role that each section plays in the proof of the theorems. This 
section is divided into two parts, corresponding to the case of the level sets of a 
single function and to the case of joint level sets. 
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2.1. Level sets of a single function. Theorem 11.11 is motivated by Question [T] 
which surfaces are homeomorphic to a component of the zero set of a harmonic 
function in M'^? Hence, to explain the gist of our approach to this theorem, we can 
restrict ourselves to this particular case. For the sake of concreteness, we shall thus 
start by trying to prove that there exists a harmonic function in M'^ whose zero set 
has a connected component diffeomorphic to the torus of genus g with N ends, as 
displayed in Figure [TJ 

The basic strategy to construct a harmonic function such that the above surface, 
which we shall call L, is diffeomorphic to a connected component of its zero set is 
the following. We will start with a function v that is harmonic in a neighborhood 
of the surface and has a level set diffeomorphic to L. We will then approximate 
this local solution by a global harmonic function u. Hence, the key for the success 
of this strategy is to be able to ensure that the level set of the local harmonic 
function is 'robust' in the sense that it is preserved, up to a diffeomorphism, under 
the perturbation corresponding to the above global approximation. 

Let us elaborate on this point to show the kind of difficulties that arise when 
carrying out this program. We have seen that the starting point of the strategy 
must be a topological stability theorem for the level sets of the local harmonic 
function v. If these level sets were compact, this could be accomplished using 
Thorn's celebrated isotopy theorem ^ Theorem 20.2]; however, the fact that the 
level sets of a harmonic function are all noncompact makes the problem much 
subtler. Moreover, this stability result must be finely tailored to provide sufficient 
control of the deformation at infinity, for otherwise this part would be a bottleneck 
in the proof of our theorems. Therefore, one of the crucial parts of the article 
is the proof of a fine C^ noncompact stability result under small perturbations, 
which holds provided that the function v satisfies a suitable gradient condition in 
a saturated neighborhood of the surface (roughly speaking, in w~^((— e, e)) in order 
to control the zero set of v) and a C^ bound. 

In view of this result, the local harmonic function v must be constructed in such a 
way that the above stability conditions hold. These conditions rule out, for example, 
an approach based on the Cauchy-Kowalewski theorem, as it does not yield enough 
information on the domain of definition of the solution. Instead, we base our 
construction on the use of Green's functions. We take an unbounded domain Q, 
whose boundary is diffeomorphic to the surface L and consider its Dirichlet Green's 
function Gq{x, y). If y is a point of this domain, Gq{-, y) defines a (local) harmonic 
function in a half-neighborhood of the boundary, which is the zero level set of this 
function and diffeomorphic to the surface L. However, the noncompact stability 
theorem cannot be applied to Go(-, y) because this function does not satisfy neither 
the saturation nor the gradient conditions. 

In order to circumvent this difficulty, we start by deforming the surface so that 
it has a controlled geometry at infinity (i.e., in this particular case, that the N 
ends of the torus are straight cylinders, as in Figure [T]). For simplicity, we keep the 
notation L := dft for this deformation of the initial surface (this is the reason why, 
without loss of generality, we are requiring rigid ends in the definition of tentacled 
(hyper)surfaces, cf Definitions 15.11 and 16.11) . We then insert a straight half- line in 
each end of the domain and, denoting by /ii the length measure on the z-th half-line. 
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It can be shown that this defines a function harmonic in the domain minus the 
half-hnes whose zero set is the 'tentacled' torus L and which does satisfy the bound- 
edness, saturation and gradient conditions of the C^ stabihty theorem in a half- 
neighborhood of the surface. Roughly speaking, the proof of this fact is based on 
estimates that exploit the asymptotic Euclidean symmetries of the construction and 
the exponential decay of the Green's function of the domain. As an aside, notice 
that to resort to this construction we will need to prove the noncompact stability 
theorem also for functions defined only in a half-neighborhood of the surface under 
consideration. 

The last ingredient is the theorem that allows to approximate the local harmonic 
function w by a global harmonic function in a certain neighborhood of the surface L. 
This approximation must be small in the C^ norm in order to apply the stability 
theorem. If L were compact, the Lax-Malgrange theorem for elliptic PDEs [501 
Theorem 3.10.7] would provide an adequate global approximation result. As L is 
noncompact, however, the situation is considerably more involved. The proof of 
our global approximation theorem relies on an iterative procedure (which does not 
apply to arbitrary elliptic PDEs) that is built over an appropriate exhaustion by 
compact sets and combines the Lax-Malgrange theorem, suitable Green's function 
estimates and a balayage-of-poles argument. 

The implementation of this strategy will be carried out in four stages: 

Stage 1: Noncompact stability theorem. In Section |3] we prove a stability theorem 
for noncompact level sets (Theorem 13.11) and a variant for functions defined only 
in a half- neighborhood of the level set (Corollarv l3.3|) . Level sets of vector- valued 
function are also treated in view of their applications to the case of joint level sets, 
to be discussed below. When applied to compact submanifolds, this result fully 
recovers Thom's isotopy theorem. 

Stage 2: Global approximation theorem. In Section |4] we prove that a local solution 
of the equation (A — q)v = 0, with q a real analytic function satisfying certain 
hypotheses, can be approximated in the C^ norm by a global solution of the equa- 
tion (Theorem 14. 5p . This extends results of Gauthier, Goldstein and Ow [12] and 
Bagby J3], dealing with equations that obey a minimum principle but not a maxi- 
mum principle. 

Stage 3: Noncompact level sets of finite topological type. In Section [5] we obtain a 
realization theorem for noncompact hypersurfaces of finite topological type, that is, 
whose homotopy groups are all finitely generated. We start by presenting the def- 
inition of tentacled hypersurfaces of finite type (Definition 15. ip and characterizing 
the hypersurfaces that can be realized as a tentacled hypersurface modulo diffeo- 
morphism fProposition 15.3]) . After a series of intermediate lemmas we show that 
for any tentacled hypersurface there exists a solution of the equation (A — k^)u = 
in M" which has a level set diffeomorphic to the hypersurface (Theorem l5.9p . Corol- 
lary \1.2\ follows from the latter theorem and Example 15.41 

Stage 4- Noncompact level sets of infinite topological type. In Section [6] we tackle 
the case of noncompact hypersurfaces of infinite topological type (Theorem l6.7p . To 
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do SO, we introduce the notion of periodic tentacled hypersurface (Definition 16. 1|) . 
discussing the associated discrete Euchdean symmetries that, together with the 
asymptotic Euchdean symmetries discussed above, play a key role in the proof 
of the theorem. Theorem IJ.il then follows from Theorems 15.91 and 16.71 a further 
generalization is presented in Remark 16. 81 

Since we have also considered level sets of the equation (A — k'^)u — 0, which can 
be compact for fc 7^ 0, for completeness it is worth discussing the case of compact 
level sets. This is done in Appendix [Aj where we show that for any compact 
hypersurface there exists a global solution having a level set diffeomorphic to it 
fTheorem lA.ip . The proof is considerably less involved as one does not have to 
deal with the lack of compactness present in Theorems 11.11 and II. 4[ and actually 
applies as well to more general equations. 

2.2. Joint level sets. We shall next deal with joint level sets of m ^ 2 harmonic 
functions, that is, with the intersection of their level sets. As we did in the case of a 
single level set, we will be interested in the case of regular level sets, which amounts 
to requiring that this intersection be regular. This leads to the consideration of 
transverse intersections of level sets uj~ (ci) n ■ • ■ fl u~^(cm), which means that the 
gradients of ui, . . . , Um are linearly independent at each point of the joint level set. 
It is well known that the topological obstruction on the submanifolds that can be 
the transverse intersection of m functions is that they must have codimension m 
and trivial normal bundle. 

Let us fist consider the case where the joint level set is compact, or, more gen- 
erally, a locally finite disjoint union of compact submanifolds of codimension m 
(with trivial normal bundle, by the above argument). The strategy we will follow 
to construct m harmonic functions with a joint level set diffeomorphic to a given 
compact submanifold L is similar to the one used for a single level set: we first con- 
struct m local harmonic functions having a joint level set diffeomorphic to L and 
then approximate them by global harmonic functions, using a topological stability 
theorem to guarantee that they have a joint level set that is also diffeomorphic to 
L. 

To construct the local solutions, we make use of the triviality of the normal 
bundle of L to characterize the submanifold as the transverse intersection of m 
hypersurfaces, which can be chosen real analytic by perturbing L a little if necessary. 
An easy application of the Cauchy-Kowalewski theorem then gives the desired 
local harmonic functions. The global approximation theorem we shall use is a 
better-than-uniform approximation result that we prove using an iterative scheme 
based on the Lax-Malgrange theorem. The key point that allows to simplify the 
treatment (and to consider more general elliptic PDEs, as in the statement of 
Theorem ll.3p is that the compactness of each component of L makes the verification 
of the topological stability conditions much less subtle than in the noncompact case; 
in particular, this is the reason why the above Cauchy-Kowalewski argument is 
successful and topological stability can been tackled using only the classical Thorn 
isotopy theorem. 

The case of joint level sets with noncompact components is more involved and, 
as in the case of a single level set, requires some fine control on the geometry of 
the submanifold at infinity. Therefore, we introduce the notion of a codimension-?7i 
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tentacled submanifold of M", which is the proper analog of a tentacled hypersurface 
and is indeed given by a transverse intersection of m such hypersurfaces. This 
automaticahy grants that the normal bundle of the submanifold is trivial. 

The strategy now is to write the codimension-m submanifold as the intersection 
of m tentacled hypersurfaces of M". One can then apply the same reasoning we used 
in the proof of Theorem 1 1.1 1 to each of the aforementioned tentacled hypersurfaces, 
obtaining m local solutions of the equations that have a level set diffeomorphic to 
one of the above tentacled hypersurfaces and satisfy the conditions of the noncom- 
pact stability theorem. Suitable control of the second-order derivatives of these 
local solutions allow us to proceed as in the case of a single level set using a global 
approximation theorem. 

Our treatment of joint level sets, which we present in Section [71 will therefore 
consist of two parts: 



Part 1: Compact joint level sets. In Subsection 17.11 we will prove that, given any 
locally finite disjoint union L of compact submanifolds in R" of codimension m and 
trivial normal bundle and a collection of real analytic elliptic differential operators 
of second order T^, there are m solutions of the equations TrUr — in M" such that 
L is diffeomorphic to a joint level set wj~^(ci) n • • • fl u~^(c„i), thereby establishing 
Theorem \1.3\ A key step in the proof is a better-than-uniform approximation 
theorem for these differential operators (Lemma 17. 2p . 

Part 2: Noncompact joint level sets. In Subsection 17.21 we will define the notion of 
a codimension-m tentacled submanifold (Definition [7?5]) and tackle the case of joint 
level sets with noncompact components, proving Theorem [77 



Notation. To conclude this section, let us present some notation that will be 
employed throughout the article. We will use the notation 

WfWcpiu) ■= max sup \D°'f{x)\ 

\a\'^PxGU 

for the CP norm of a map / in a set t/ C M", and denote by B{x,p) the ball of 
radius p centered at a point x £ R". We will always assume that n ^ 3. If 51 is a 
(possibly unbounded) domain of R", we will denote by 



inf I f |Vwp dx:ve C^i^), f v^ dx = l\ 



the infimum of the L^ spectrum of —A in the domain fl with Dirichlet boundary 
conditions. 

As customary, we will say that a function satisfies a PDE in a closed set S if the 
PDE holds in some open set containing 5*. All the diffeomorphisms appearing in 
this article are assumed C°° and orientation-preserving without further mention, 
and all the submanifolds of R" that we will consider are C°° , oriented, properly 
embedded and without boundary. 
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3. NONCOMPACT STABILITY THEOREM 

In forthcoming sections, it will be essential to ensure the stability under small 
perturbations of certain level sets of various maps. When these level sets are com- 
pact and regular, this is granted by Thom's isotopy theorem [U Theorem 20.2], but 
noncompact regular level sets are generally not stable. 

Our goal in this section is to prove the following C^ stability theorem for non- 
compact level sets, which is of separate interest and allows us to deal with uniform 
perturbations of the type considered in Section U) Our proof is based on totally 
different ideas from those of Thom and has the crucial advantage of being explicit 
enough to allow for a fine control of the diffeomorphism at infinity. In our proof, the 
diffeomorphism is constructed essentially as the time-1 flow of a carefully chosen 
vector field, whose components are checked to be suitably small: 

Theorem 3.1. Let U be a domain in K" and let f : U ^ M™ be a C°° map, with 
m < n. Consider a (possibly unbounded) connected component L of the zero set 
/~^(0) and suppose that: 

(i) There exists a domain V whose closure is contained in U and such that the 
component of f~^(B{0,ri)) connected with L is contained in V for some 
77 > 0. 
(ii) The gradients of the components fr of the map f satisfy the condition 

inf |wiV/i(a::)-f •••+w„V/™(a:)| >0 

for all vectors uj G R™ of unit norm. 

Then, given any e > 0, there exists some 5 > Q such that for any smooth function 
g-.U-)- M'" with 

(3.1) \\f-9\\cHv)<S 

one can transform, L by a diffeomorphism, <I> o/M" so that $(i) is the intersection 
of the zero set g~^{0) with V . The diffeomorphism $ only differs from the identity 
in a proper subset of V and satisfies ||$ — idllfjomn-) < e. 

Moreover, if the C^ norm \\f — g\\cp(v) ** small enough and the first p + 1 
derivatives of the maps f and g are bounded (i.e., \\f\\cp+'^(V) O'f^d \\g\\cv+'^(v) O''"'^ 
finite), then the diffeomorphism is C^-close to the identity: ||$ — id||cp(Rii) < e. 

Before presenting the proof of the theorem, a comment on its hypotheses is in 
order. Condition (i) asserts that V contains a neighborhood of the submanifold L 
saturated by the map /. Condition (ii) imposes a gradient bound on the function / 
in the set V, and was also used by Rabier in his extension of Ehresmann's fibration 
theorem [22] (notice, however, that the ideas we use in our proof are not related 
to Rabier's). When m ^ I, this condition simply asserts that |V/| > C > in 
V, while when m > I this condition measures to what extent the gradients of the 
components of the map / = (/i, . . . , /„) remain linearly independent in the set V. 

Proof of Theorem \3.1\ By (i) we can assume that the open set V is exactly f^^{B{0, rj)) 
and C/ is a small neighborhood of the closure of V where condition (ii) and ()3.ip still 
hold. Let us take an open interval / D [0, 1] and define the map h : I xU ^^ M.x M™ 
as 

(3.2) hit,x):={t,{l-t)fix)+tgix)). 
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Thus h{t, x) is an auxiliary map which connects the functions / and g in the sense 
that /i(0,x) = (0,/(a;)) and h{l,x) = (1,5(2:)). 

Let us take an arbitrary point x in U and denote by F the m x n matrix Df(x), 
with components Fij :— djfi(x). Denoting by F* the transpose matrix of _F, by 
the condition (ii) it is clear that the lower bound 

(3.3) \F*^\ = |wiV/i(a;) H h WmV/m(a;)| ^ ci\uj\ 

holds for all vectors w G M™, where ci is a positive constant independent of the 
point X G U . 

Our goal now is to show that the derivative of the map h{t,x) is uniformly 
bounded from below. Consider the {m + 1) x {n + 1) matrix H := Dh(t, x), which 
we write as H = Hq + A, with 

/10\ . _ ( 

"^"lo F)^ '^'-\g{x)^f{x) t{Dg{x)^F) 

Here we are using the notation D, V for the derivatives with respect to the variables 
t and X, as opposed to the derivatives i^,V with respect to x. Setting ||yl|| := 
sup|„|^i \Av\, it is clear that ||^|| ^ C5 by the C^ estimate (|3.1I) . Since 

\H*v\^\\H*v\-\A*v\\ 
and \A*v\ ^ ||j4|||u| ^ C(5|u|, it stems from the inequality p.3p that 

(3.4) \H*v\ ^ (min{l, ci} - C(5) \v\ > C2|t;| 

for small enough 5 and all w G R x K™, (i, x) € I xU. Here C2 is a positive constant 
independent of t and x. Therefore, the kernel of H* is trivial, so H is onto and 
thus the map h (resp. /) is submersive in I x U (rcsp. U). 

Let us denote the components of the map h as (ho, hi, ... , h„i), with /io(i, x) = t 
and hi{t,x) — (1 — t) fi{x) +tgi{x) for 1 ^ i ^ to. The condition p.4p ensures 
that the self-adjoint (to, + 1) x (to, + 1) matrix HH* satisfies HH* ^ c^, so that its 
inverse B :— {HH*)~^ is positive definite and satisfies the inequality B ^ C2^ for 
all (i,x) £ I X U. We will denote by {B^v)^^,^^ the matrix elements of B. 

In order to construct the desired diffeomorphism, for each integer /j, = 0, 1, . . . , to, 
it is convenient to consider the vector field S^ defined by 

rri 

(3.5) 5^:=^B^AV/iA(t,a;) 

A=0 

in / X [/. By definition, the scalar product (in M"+^) of these vector fields with the 
space-time gradient of the function h^ is given by 

ni m n 

(3.6) Sf, -VK^Y. ^mA V/iA ■^K^Y.Y. B^.xHxpH,p = 6^, , 

A=0 A=0^=0 

with 6f^i, standing for Kronecker's delta. Hence these vector fields define a paral- 
lelism; in particular, each field S*^ does not vanish in / x [/. Besides, denoting by 
(SfiO, Sfj,i, . . . , S^n) the components of S*^, one can readily check that the compo- 
nent Sfj,i3 coincides with the matrix element (_BiJ)^^, so that the norm of these 
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vector fields is bounded in I x U by 

n 

(3.7) |5^p = 5]52^ = i?^^^c-^ 

/3=0 

Choosing a small enough positive constant (5 < ^ and using the condition (|3.ip , 
we can ensure that g^^{0) is contained in the proper subset W :— f^^{B{0, ^)) of 
U. Let us take a C°° function x ■ ^ >^ U ^- [0,1] which is equal to 1 in [0, 1] x 14^ 
and is supported in I x U. 

From Eq. (I3.6p with indices ^ — ly — it stems that one can decompose the 
vector field 5*0 as 

(3.8) SQ=:dt+X, 

where X is a vector field in I x U orthogonal to the coordinate vector field dt ■ Let 
us now consider the vector field 5*0 := dt + X-^i whose time-1 fiow will essentially 
yield the desired diffeomorphism. Since l^ol < l^o] is bounded in I x U by the 
estimate l\3.7\i and coincides with the coordinate field dt in a neighborhood of the 
boundary d{I x U), it follows that the vector field Sq defines a C°° local fiow (ps 
in I X U. Let us write this local fiow as 

(3.9) ip,{t,x)=:{t + .3,^s,t,x)) 

and notice that, for any x ^ U, (ps{0, x) is well defined at least for s £ [0, 1]. 

Let us define the codimension-(TO + 1) submanifolds 

Lo := {0} X f-\0) , L, := {1} X g-\0) 

of the set I xW. Clearly Lq = h-^{0,0) and Li = h-^{l,0). We claim that the 
time-s fiow (ps{Lo) of Lq is contained in the time slice {s} x W for all s G [0, 1]. In 
order to see this, let us take (0, xq) G Lq and use the notation (s, x(s)) :— (/3s(0, a;o). 
As So ■ S7hi = for 1 < j < m by Eq. ([SJ]) . one finds that 



— h^{s,x{s)) 



il~x)\X-Vh,is,x{s))\^il-x) 
{I - x)\Mx{s)) ~ 9^{x{s))\ < 6 



^(^,-(^)) 



by the estimate (|3.ip . so that \hi{s, x{s))\ < S for all s G [0, 1]. In turn, this yields 
the bound 

\Mx{s))\<2S 

for all s G [0; 1] and 1 ^ i ^ ?7i, so that the time-s fiow ipsiLo) is contained in 
[0, 1] X VK for s G [0, 1]. Since x = 1 in [0, 1] x W, actually h,{s,x{s)) = h,{0,xo) = 
0, which implies that (pi(Lo) C Li. By reversing the argument one infers that 
(p_i(Li) C Lq, so that (pi(Lo) = Li, as claimed. 

Let us now consider the map $1(2:) :— $(1,0, a;) (see Eq. (|3.9p l. which is a 
diffeomorphism of U by construction. As the map $1 : U -i- U so defined is then 
equal to the identity in a neighborhood of dU, it can be trivially extended to a C°° 
diffeomorphism $ of M" by setting 

$(.):= |*^(^)' ^^^' 
^ ^ \x, x(^U. 

By construction, $(£) coincides with g^^(O) n V. 
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To complete the proof of the theorem, it only remains to show that the diffeo- 
morphism $ is close to the identity. Since |x| ^ 1 and $i is defined in terms of the 
local flow (I3.9p . in order to show that the norm |j<i> — id||(^o(Rii) is small it suffices 
to prove that ||-'^||c''(/x(7) is close to zero, where X is the vector field defined in 
Eq. ()3.8|) . Since the space-time gradient of /iq is 9t, it follows from the definition of 
X and Eq. (|3.6p that X can be written as a linear combination of the gradients of 

X = ^v^,Vh^. 

11=0 

Here the function w^ = f^(i, x) depends on the variables t and x. Notice, moreover, 
that 



X-Vh^^So- Vhf, 



by Eq. (I3.6P . Hence the norm of the vector field X satisfies 



II ^ I-^-EIILo^m^^mI EI"i h\ l/» -g'l < TO|i^|sup^|/-g| rnS^ 

by Eqs. p.ip and p.4p . thereby proving the smallness of the C° norm ||$— id||(70(K>i). 

To estimate the C^ norm of $ — id, one argues as in the case of the C° norm but 
one needs to estimate the C^ norm of the vector field xX. First of all, it should be 
noticed that the hypothesis that the derivatives of / are bounded ensures that the 
set V contains a 'metric tubular neighborhood' 

{x e R" : dist(a;, L) < p} 

of the submanifold L, where p is some positive constant. This obviously allows us 
to assume that all the derivatives of the function x a-re bounded in M". After a 
straightforward but tedious computation, one finds that the C^ norm of the vector 
field X can be controlled as 

\DPX\ ^ C7(||/||c.+i, ll.9llcp+i, llxllc) 11/ - 9\\c^ 

in the set I xU, where Cdl/Hc-p+i, ||ff||cp+ii llxllcp) is a constant depending on the 
norms oi f, g and x- The theorem then follows. D 

Remark 3.2. When the submanifold L is compact, all the derivatives of the maps 
/ and g are necessarily finite in a suitable neighborhood of L, so the C condition 
11/ ^ 5llcp(V) < S automatically yields a diffeomorphism $ that is close to the 
identity in the C^ norm (as was proved in Thom's isotopy theorem). 

As a matter of fact, it is worth pointing out that the above proof does not only 
demonstrate Theorem l3.11 but also the following closely related result, in which the 
component L of the zero set of the function is allowed to lie on the boundary of the 
set V. Although this is only a minor modification, we will find this result of use in 
forthcoming sections: 

Corollary 3.3. Let U be a domain in M" and let f : U ^- K™ be a C°° map, 
with m < n. Consider a (possibly unbounded) connected component L o//^^(0). 
Suppose that: 
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(i) There exists a domain V whose closure is contained in U and such that 

the intersection /~^(0) ClV is L. 
(ii) The gradients of the components fr of the map f satisfy the condition 

inf IwiV/i + • ■ • + w„, V/„j| > 
for all unit vectors uj G M™ . 

Then, given any e > there exists some 5 > such that for any sm,ooth function 
g : U ^- K™ with \\f ~ g\\c^{V) < ^ one can transform, L by a diffeomorphism $ 
o/K" so that $(i) is the intersection of the zero set g~^{Q) with the closure of V 
provided that any component of the set 

(3.10) {x<eU ■.{l-t)f{x)+tg{x)^Q} 

that intersects V is contained in V for all t G [0, 1]. Besides $ — id is supported in 
a proper subset of any fixed neighborhood ofV and \\<^ — \d\\co(m.n) < e. If moreover 
the C'P norm ||/ — g\\cp(v) *'^ small enough and the first p + 1 derivatives of the 
maps f and g are bounded in V , then one can take the diffeomorphism $ C^-close 
to the identity. 

Remark 3.4. In this paper we will only apply this corollary to g :— f ~ c, in which 
case the sets p.lOp are simply f~^{tc). Clearly a straightforward modification of 
Eq. p.2p allows to deal with the case where the level sets f~^{tc) do not have the 
above property but there is a smooth curve 7 : [0,1] -^ B{0,S) with 7(0) — 0, 
7(1) = c and such that any component of f~^{j{t)) meeting V is contained in V. 
An analogous result can also be established for the case of two arbitrary maps / 
and g. 

4. A GLOBAL APPROXIMATION THEOREM 

Our goal in this section is to establish that any function v satisfying {A—k'^)v = 
in a closed (possibly unbounded) set S C R" can be approximated in the C^ norm 
in this set by a global solution of the latter equation, provided that the complement 
of the set S does not have any bounded connected components. This condition on 
S is used to apply the Lax-Malgrange theorem and in an argument on the balayage 
of poles. 

We will prove this result for the equation 

iA-q)v^O in 5 C R" , 

as the proof for {A — k'^)u — is not substantially easier. Here g is a bounded, non- 
negative function on R" which we take real analytic (this condition is not necessary, 
as all the arguments in this section can be modified to deal with lower regularity, 
but this allows us to avoid inessential technicalities). In order to get a C^ uniform 
approximation result, we will assume that the C^~^ norm H^Ucp-^m") is finite. 
The results we prove in this section generalize theorems of Gauthier, Goldstein and 
Ow '12] and Bagby |3j to equations that satisfy a minimum principle but not a 
maximum principle. 

Throughout the paper, we shall often need the Dirichlet Green's function Gn {x, y) 
of the operator A — g in an unbounded domain 17. For completeness, in the following 
proposition we prove the existence of a minimal Dirichlet Green's function using 
the classical method of compact exhaustions and summarize some properties that 
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will be required later. When fi = M", we will write G{x, y) instead of Gk" [x, y) for 
the ease of notation. 

Proposition 4.1. Let Q be a domain in M", possibly unbounded and with nonempty 
boundary of class C°° . Then there exists a minimal positive Dirichlet Green's func- 
tion 

Gn'-i^lx f])\diag(f7 xfl) ^R, 

which satisfies: 

(i) The Green's function is symmetric, that is, Gsi(a;,y) — Gi-i{y,x). 
(ii) Gn(-,y) G C°° {Ti\{y}) r\ Ll^{n) and {A^ - q{x))Gn{x,y) = -5y{x), with 

6y being the Dirac measure supported at y. 
(iii) Gsi(-, y) = on the boundary Qfi and Gn{x, y) ^ G\x — y\'^^"', with G^^ = 
(n-2)|§"-i|. 

Proof. If il is bounded, the result is well known, so we henceforth assume that fl 
is unbounded. Let ili C J^i C 1^2 C ■ • • C fi be an exhaustion of Q by bounded 
domains with smooth boundaries and let 

G, : {^i X f^j)\diag(rJi x n,) -> R 

be the Dirichlet Green's function of the domain fli, which is symmetric and satisfies 

(Aa; -~q{x))Gi{x,y) = -5y{x) 

in Vti and Gi{-,y)\an, == 0. 

Since g ^ 0, the minimum principle ensures that the Green's function Gi{-,y) is 
positive in Vli and that the sequence of Green's functions is monotonically increasing 
in the sense that 

(4.1) G^{x,y) <G^+i{x,y) 

for all x,y G fli- Moreover, the minimum principle also guarantees that Gi is 
bounded by the Green's function of the Laplacian in W as 

G,{x,y)^C\x-y\'-^, 

with G^^ := {n — 2)|S"^^|. By the monotonicity property (|4.ip . (^Gi{x,y)) , is 
then a Cauchy sequence. Standard gradient estimates [13j imply that, for any open 
sets B CB C B' cW c{n,r\ nj)\{y}, 

||G.(.,y) - G,(-,2/)||c^(s) < G'||G.(-,y) - G,(-,y)||co(B') , 

where G' only depends on B, B' and on the function q. Hence Gi[-,y) converges 
G^-uniformly on compact subsets of n\{y} to a solution Gni-^y) of the equation 
(A — q)Gn{-,y) — ^Sy. The fact that Gn has the properties (i)-(iii) stems from 
the above construction. D 

The proof of the approximation theorem relies on the following three lemmas, 
whose proofs make use of several ideas of Bagby [3] and Bagby and Gauthier [3] 
but rely on a different argument because the differential equation (A — q)v = 
does not satisfy a maximum principle. Lemma 14.21 estimates the G^ norm of the 
difference between Dirichlet Green's functions with nearby poles using Schauder 
estimates (throughout this section, p will denote an arbitrary positive integer): 
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Lemma 4.2. Let V be an open subset o/M" and let y be a point in V. Then for 
any e > there exists an open neighborhood By d V of y such that 

\\G{-,y)-G{-,z)\\cP(R^\v) <e 
for all z G By . 

Proof. We can assume that V is bounded without loss of gencrahty. Let us arbi- 
trarily fix some real constant p > 0. Since (A — q)G{-,z) — in M"\{z}, interior 
Schauder estimates [TS] Corollary 6.3] show that there exists a constant Ci (de- 
pending solely on n, p, p and ||'7||cp-2(]jn)) such that the C^ pointwise estimate 

(4.2) max {D^Gix, z)\ ^ Ci||G(-, z)||co(b(.,p)) < G2{\x - z\ - p)'~" 

holds whenever jx — z| > p. Here the last inequality follows from the third statement 
in Proposition [131 thereby showing that D°'G{-,z) tends to zero at infinity for all 
\a\ <p. 

By the G^ estimate (|4.2p and the boundedness of the set V, for any e > one 
can take a compact subset K D V oi M" such that 

(4.3) sup sup \D'^G{x,z)\ < - 

for I a I ^ p. Moreover, since D^G{x,z) depends continuously on z G R"\{a;} and 
K is compact, there is a small neighborhood By C V oi y such that 

sup sup \D"G{x,y) ~ D"G{x,z)\ < e. 
zeBy xeK\v 

By the definition of the set K, the lemma follows. D 

In the following lemma we show that the Green's function with pole y can be 
approximated in a suitable sense by finite linear combinations of Green's functions 
with poles in a prescribed set. The proof is based on density and duality arguments: 

Lemma 4.3. Given an open set W C M", let us take a compact subset K C W of 
nonempty interior and a point y G W . Then, for any e > 0, there exist a finite set 
V <Z K and constants {az)zev C M such that 






< e 

Cp(R"\W) 



Proof. We assume that the point y does not belong to the set K, since otherwise the 
statement is trivial, and consider a proper open subset W C W containing y and 
-ftT, with smooth boundary and such that p '■= ^ Aisi{dW, dW) > 0. (Later on, we 
will take advantage of this set to transform a C° estimate into a C' estimate.) Let 
Gq{M."\W') denote the Banach space of continuous functions on the complement 
of W tending to at infinity, endowed with the supremum norm. We denote by 
K, the subspace of Go{W-\W') consisting of all finite linear combinations of G(-, z) 
with z € K. 

By the generalized Riesz-Markov theorem, the dual of Co(M"\W^') is the space 
M(M."'\W') of the finite signed Borel measures on K" whose support is contained 
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in the complement of W' . Let us take any measm^e /i G J^(W"-\W') such that 
J f dfi = for ah f £ K.. Let us now define a function F g Ll^^{R^) by 

F{x) ■- I G(z,x)dn{z), 

so that F satisfies the equation (A — q)F = —fi. Since F is identically zero on the 
set K by the definition of /i and K has nonempty interior, the unique continuation 
theorem ensures that the function F vanishes on W' . It then follows that ^ also 
annihilates G{-,y) because 



F{v)^ j G{z,v)dii{z) 



-0, 



which shows that G(-, y) can be uniformly approximated on IR"\VF' by elements of 
the subspace /C as a consequence of the Hahn-Banach theorem. 

To complete the proof of the theorem, let us take a sequence (/i)^i C K, such 
that 

\\G{-,y) -/j||co(K"\H") < -, 

which is guaranteed to exist by the above argument. We can now use the Schauder 
estimate (14.21) to show that, for all x G M"\Vl^, one has 

max|D^G(x,y)-D"/,(a:)|^Ci||G(-,y)-/,||co(B(,.p))< — 

with Ci independent of x, as claimed. D 

In the following key lemma we discuss the balayage of the poles, which basically 
allows us to get rid of the poles by taking them to infinity, exploiting the fact that 
the set V is unbounded. The proof relies on an iterative argument that utilizes 
Lemma l473l to sweep the poles further in each step. 

Lemma 4.4. Let V he an unbounded domain of W^ containing a point y. Then, 
for any e > 0, there is a function w which satisfies the equation (A — q)w — Q in 
R" and approximates G{-,y) as 

\\G{-,y) -w||cp(R"\y) < e. 

Proof. As the set V is unbounded, we can take a parametrized curve 7 : [0, 00) -^ V 
without self-intersections (that is, 7 is injective) such that 

7(0) = y, lim |7(i)| = 00. 

r— ^00 

This is the curve that we will use to 'sweep the pole to infinity'. For each nonneg- 

ative integer i we denote by Ki C V Ci B{'-f{i), ^) a compact neighborhood of the 

point 7(i), and we let {Wi)°^i C F be a family of bounded domains such that each 

Wi contains both A'j_i and Ki. We also require that Wi is 'narrow enough' in the 

sense that 

sup dist (x, 7([i — 1, i])) < 1 . 
xeWi 

Moreover, we introduce the notation /Cj for the space of finite linear combinations 

of G(-, z) with poles z E Ki. 

Let us fix e > 0. By Lemma [4.31 there exists wi G /Ci such that 
l|G(-,y) -wi||cp(R"\Wi) < o ■ 
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Since wi is a finite linear combination of Green's functions with poles in Ki , using 
Lemma 14.31 one can inductively show that there are functions Wi £ ICi such that 

(4-4) \\wi - Wi_i||cp(R"\vi/.) < ^ , 

for all i ^ 2. Since the distance between the point y and the set Wi tends to 
infinity as i — > oo, this shows that Wi converges C^-uniformly on compact sets to a 
function w that solves the equation (A — q)w = in R". The lemma now follows 
upon noticing that 

i 

\\G{;y)-m\\c,^K^\v) ^ l|G(-'y)~^i|lcp(R'.\v)+I^II"'j'""'j-i|lcp(R"\n 

i 

for all j, on account of the C^ estimate (14.41). D 



Armed with the previous lemmas, we are now ready to prove the main result of 
this section, which will be frequently used in the rest of the paper. The rough idea 
of the proof is to extend the local solution w to a smooth function w in M" supported 
in a neighborhood of 5"; using an iterative procedure based on Lemmas 14. 31 and 14. 41 
vj is then approximated by a linear combination of Green's functions with poles in 
the complement of S and subsequently the poles are swept off to infinity. 

Theorem 4.5. Let S be a closed subset of M" whose complement does not have 
any relatively compact connected components. Then any function v that satisfies 
the equation (A — q)v — in S can be approximated in the C^ norm by a global 
solution to this equation. (That is to say, for any e > there exists a function u 
satisfying (A — q)u = in M." with \\u — v\\cp{S) < ^-J 

Proof. By hypothesis, there exists an open subset fl D S such that (A — q)v = 
in n. Let us take a smooth function x ■ I^" ^ ^ equal to 1 in a neighborhood of S 
and identically zero outside fi, and define a smooth extension of v to M" by setting 
w := xu. 

Consider an exhaustion of K" by bounded domains fli as in the proof of Propo- 
sition HTTl which we choose so that M"\rii does not have any compact components. 
We define the sets Ui := 17i_|_i\rii_i, with J7o := 0, and take smooth functions (pi on 
M" whose support is contained in Ui and such that X]i=i "ySi = 1 in a neighborhood 
of Qj . A possible way to construct these functions is to take smooth functions 0i 
compactly supported in fli-^-i and such that (pi = 1 in a neighborhood of fli] one 
can now set ipi := (pi — (pi-i, with (f>o '.= 0. 

We will now write the function w as a sum of global solutions to the equation 
and of functions that satisfy the equation but in a compact set. To do so, for each 
positive integer i let us define the smooth function 

f, := Lpiiq- A)w, 

which is supported in {^\S) D Ui by the definition of the functions x ^nd ipi. 
Defining the function 

Wi{x) := / G{x,y) f,{y)dy, 
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it is apparent that Wi satisfies the equation (A — q)wi — —fi, so that 

in ^j. By applying the Lax-Malgrange theorem (201 Theorem 3.10.7] inductively, 
one easily infers that there exist functions gi satisfying the equation (A — q)gi ~ 
in ]R" and such that 



(4.5) 



w - 



1 
< - 



for any positive integer ]. Therefore the function w can be expressed as 

oo 

(4.6) w = ^(w, +3,), 

the convergence being C^-uniform on compact subsets of R". 

Let us now approximate Wi by functions wi that satisfy the equation in M" minus 
a finite set of points. Let us use the notation ||/i||Li := / l/j(2;)l 'i-^ for the L^ norm 
of ]i and suppose that this quantity is nonzero. By Lemma 14. 2[ for each point y 
in the support of the function /^ one can find a neighborhood By 3 y contained in 
(M"\S') n U, such that 

\\G{-,y) - G{-,z)\\cp{su(R"\Ui)) < 11. II ^2i+i 

for all z G By. As the support of fi is compact, there exist finite pairwise disjoint 
sets Vi C supp fi such that supp fi is contained in the union of balls Ui/eP ^v ^^^ 
{Xv)v&Vi C C^{M.^) be a partition of unity subordinated to the balls {By : y G Vi}, 
that is, nonnegative functions such that supp Xy C By and 

for all X in a neighborhood of supp fi. Defining the function Wi as 






w,\x) := 2^ Lr(x, J/) / y„\z] Ji[z)az 
ydVi 

and using that Xy i^) defines a partition of unity, one easily obtains the C^ estimate 
\D^w,{x) - D"w,{x)\ = Y. j {DZG{x, z) - DZG{x, y)) Xy{z) f^[z) dz 



(4.7) < 



yeV: 
e 



T^m y l/»WI ( Y. ^!^W) ^^ = ^ 



\\f^\W- . ^y^-p^ 



which holds pointwise in S* U (M"\C/i) for all |a| ^ p. When ||/i||Li = 0, we can 
simply set Wi := 0. 

We will next use a balayage argument to sweep the poles of Wi. As R"\S' does 
not have any bounded connected components, one can take (possibly disconnected) 
pairwise disjoint, open unbounded sets Vi C IR"\S' such that each point y € Vi is 
contained in an unbounded component of Vi . Notice that the sets Vi can be chosen 
so that for each compact set K CM.'^ there exists an integer J such that K CiVi = 
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for all z ^ J. As Wi is a finite linear combination of Green's functions with poles in 
supp/i C Ui\S, from Lemma l4.4l it follows that there exists a function hi satisfying 
the equation (A — q)hi = in M" and such that 



(4.8) 



\Wi - /li||cp(R"\y,) < 



2^+1 • 
As we will see, this condition ensures that the function 



(4.9) 



" ■='^i'^t+9i) 



i=l 



is well defined and that the sum converges C^-uniformly on compact subsets of 
M". In order to show this, let us take an arbitrary compact subset K of R" and an 
integer J such that ViCiK = UiDK — ^ ioT all i ^ J , which is known to exist by the 
way the sets Vi and Ui have been defined. In particular, the latter condition implies 
that K is contained in ilj-i- The C^-uniform convergence of the series (|4.9p on 
the set K then follows from the pointwise C^ estimate on K 



i=j 



Y^ {D^h, + D^'g,) ^ Y^ \D"hi - D°'w^\ + ^ \D°'w^ - D°'wi\ 

i—J i=^J 

k 

k 

i=l 
.1-1 

Y {D^gi + D'^w 



^ \D°'w\ 



J-1 



J2 {D''g^ + D^W,) 



i=l 



+ e, 



s$ iD^wl 



k 



where \a\ ^ p, k > J and we have used Eqs. (|4.5p - (|4.8p . As a consequence of the C^- 
uniform convergence on compact sets, clearly u satisfies the equation (A — q)u = 

in R". 

Finally, one can invoke the uniform convergence of the series (|4.6p and (|4.9p to 
show that u approximates v in the set S in the C^ norm: 



OO 

E 



"i\\CP{S) 



i=l 



OO 



z=l 



4=1 j=l 



as we wanted to prove. 



D 



5. NONCOMPACT LEVEL SETS: THE CASE OF FINITE TOPOLOGICAL TYPE 

Our goal in the following two sections is to study which hypersurfaces can be 
level sets of solutions to the (Laplace or Yukawa) equation 

(A - k^)u = 
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in K", with k a real constant. For convenience, in this section we wiU focus on the 
case where the hypersurfaces have finite topological type, postponing the case of 
infinitely generated hypersurfaces until Section [51 

In this section we will establish a sufficient condition that allows us to show that 
a wide class of hypersurfaces can be realized as level sets of a solution to the above 
equation modulo diffeomorphism. As mentioned in Section [21 a key step in the 
proof is to exploit the 'freedom' associated to this diffeomorphism to embed the 
hypersurface in a way that permits us to have a fine control at infinity of various 
quantities. These suitably embedded hypersurfaces will be called tentacled. 

It is standard that the notion of infinity in a hypersurface is captured in a 
precise way by its end structure. (Let us recall that, roughly speaking, an end of 
a noncompact manifold L is a component of L\K for a sufficiently large compact 
subset K C L; for the precise definition, cf. \W). An end is said to be (smoothly) 
collared if it has a neighborhood diffeomorphic to S x [0, oo), where E is a compact 
submanifold of codimension 1 in L. Tentacled hypersurfaces (of finite type) are 
simply hypersurfaces with a finite number of ends, which are all collared and whose 
geometry is suitably controlled: 

Definition 5.1. An unbounded domain fi of M" with smooth boundary of finite 
type is tentacled if one can find J embedded images of K"^^ in M", which we will 
call IIj, such that the following statements hold (cf. Figure [5]): 

(i) Each Hj divides R" into two domains, H^ and HJ^ , and the sets H^ , • . ■ , -ff j 
are pairwise disjoint. 

(ii) The intersection J7 n H^ n • ■ • fl HJ is bounded. 

(iii) Each connected component of the (possibly disconnected) hypersurface 
with boundary il n (Hi U • • • U IIj) is bounded and contained in an affine 
hyperplane. These components will be denoted by A^, with 1 ^ i ^ N. 
Here N denotes the number of ends of the domain. 

(iv) Each 'tentacle' Si, which is the component of fir\{H^Li- ■ -LiHj) connected 
with the 'cap' A^, is isometric to the Riemannian product A^ x (0, oo); in 
particular, the closure of A^ intersects orthogonally the boundary of il. 

A (possibly disconnected) hypersurface of finite type is tentacled if it has a finite 
number of connected components, each of which is the boundary of a tentacled 
domain. 

Remark 5.2. As we have allowed the intersection fi n Hj to be disconnected, the 
different tentacles Si whose caps A^ are contained in the same intersection fl n Hj 
can be linked among them. In this case, the number of ends N is larger than the 
number of hyperplanes J. 

A very wide class of hypersurfaces of R" can be transformed into a tentacled 
hypersurface via an appropriate diffeomorphism of M" . In the following proposition 
we will characterize the class of hypersurfaces that are equivalent to a tentacled 
hypersurface modulo diffeomorphism. Although the proof is elementary (and will be 
safely omitted), this result is of interest because it automatically provides a number 
of nontrivial examples of hypersurfaces diffeomorphic to a tentacled submanifold, 
which will be subsequently shown to be diffeomorphic to a level set of a harmonic 
function in M": 
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Figure 2. A tentacled surface L diffeomorphic to the torus with 
iV == 3 ends. 



Proposition 5.3. Let M he a noncompact connected hypersurface o/M". Then 
there exists a dijjeomorphism $ o/R" transforming M into a tentacled hypersurface 
$(M) if and only if there is a compact hypersurface with boundary M. of the closure 
of the n-ball B" such that: 

(i) The intersection of M. with the boundary of the n-ball is precisely its bound- 
ary dA4, and this intersection is transverse. 

(ii) There is an embedding /i : M" ^ B" such that h{M") = B" and h{M) is 
the interior of M. . 

Example 5.4. By Richards' classification theorem [23], it is straightforward that 
any noncompact surface of finite type (i.e., the genus g torus with N ends considered 
in Question [1]) can be embedded as a tentacled submanifold in E'^. More generally, 
a theorem of Calcut and King [7] ensures that any real algebraic hypersurface of M" 
satisfies the hypotheses of Proposition 15.31 and can thus be realized as a tentacled 
hypersurface modulo diffeomorphism. 

Let us begin by proving some intermediate lemmas. To construct the local 
solution of the equation (A — fc^)u = that will be subsequently approximated by 
a global one, a basic tool is the Dirichlet Green's function Gn{x, y) of the operator 
A — fc^. Therefore, we will make use of a pointwise estimate of the Green's function 
and its second-order derivatives that for our purposes is most conveniently stated 
as follows. Up to first order, this estimate is proved in Ref. [211 Theorem 5.7]; in the 
case of second-order derivatives, the argument follows the same lines and was kindly 
communicated to us by Yoichi Miyazaki (more generally, it yields C™ estimates for 
the Green's function of appropriate elliptic operators of order to). For the benefit of 
the reader, we will sketch the proof of this result below. The definition of a uniform 
C domain is given in [26. Definition 3.2]; in this paper, we will only need to know 
the (quite evident) fact that tentacled domains (or periodic tentacled domains, to 
be defined in Section |6] below) are uniform C domains for all r. 
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Theorem 5.5. Let fl be a uniform C* domain of M". Then, for any positive 
constant e we have the pointwise estimate 



\D^D^^Gnix,y)\^C\, 



y 



2-«-|Qh|/J| C'(e-(An+fc=^)i/2)|2;-y| 



for the Green's function and its derivatives, which holds for all points {x,y) E flxQ 
and any multiindices with \a\,\f3\ ^ 2. Here Xq is the bottom of the spectrum of 
—A in the domain fl with Dirichlet boundary conditions and C, C are positive 
constants. 

Proof. When \a\, |/3| ^ 1, this resuh is proved in 26j Theorem 5.7] using the key 
estimate [3S1 Theorem 5.5] 

(5.1) P"(Ao - A)-iiP^iP ^ C\\n + A|^-i , 

vahd for \a\ ^ 2. Here A is a possibly complex constant (whose real part we assmxie 
to be larger than — Ao), || • \\lp^li stands for the operator norm LP{VI) — s> L''{Vl) 
and we are denoting by Aq the Laplacian in Vi with Dirichlet boundary conditions. 
Hence in what follows we will show how the argument can be modified to allow 
multiindices with |a|, |/3| ^ 2. 

Throughout this proof, for the sake of simplicity we will denote by C generic 
positive constants independent of A. The L^ — ?• L'' norm of the resolvent can be 
readily obtained from the above inequality and the Sobolev inequality 

(5.2) Iklk. ^C||u;||i;^^^~'VllJit'\ 
thereby finding that 

||(Ao - A)-1||lp^l, < C\\n + Api+^'i-^) 

whenever 1 < p < q < oo and p~^ — q^^ < 2/n. 

It is a trivial matter to see that the L^ — > L^ bound for the resolvent (j5.1|) yields 
the inequality ||(Ao — X)^^\\lp^lp ^ CjAo + A|^^, while by elliptic regularity it is 
well known that \\{Aq — X)^^\\iyi.p_^w^-p ^ C*. Hence once obtains that 

||(Ao-A)-2||^p^^3,p <; \\{An-Xr^\\Lp^w^:p\\{^n-Xr^\\w^.p-,w^.p i^C\Xn+X\-^ 
Therefore, using these estimates and the Sobolev inequality 



\D^w\\l^ ^ C\\wh,' 



w 



Vya.p 



for |a| ^ 2 and p > n wc immediately get the L^ ^ L°° bound 
(5.3) p"(Ao - Xr'hp^L- ^ C\Xn + Ap^+ld^l+i^) 

for the above range of parameters. 
The next step is to prove that 

IKAf, - X)-^\\l2^lp < C\Xn + Ar"+t(i-i) 

for all finite p > n and large enough m. To this end, it suffices to take m sufficiently 
large, so that there are numbers 2 =: po < pi < ■ ■ ■ < pm '.— p whose consecutive 
inverses satisfy p~\ — p^^ < f • Then one can apply the L^ -^ L'^ estimate ()5.2p 
to each pair of consecutive indices to yield the desired L^ — > L^ bound. 
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Now we can combine the above equation and the estimate (|5.3p to derive that, 
with in := i ~ 2, 

||D"(Ao - \)-'\\v^-^L^ ^ C\\n + Ar-'"+^ + t 
for \a\ ^ 2. Since the image of (Ao — A)^-' is obviously contained in the Holder 
space C ^~p provided j is chosen as above, we can apply Tanabe's ST* lemma [2£| 
Lemma 5.10] with S := D°'{Aq - A)"-'', T := Df{An - X)~^ and multiindices with 
|q;|, \/3\ ^ 2 to derive that the integral kernel K{x,y) of the operator (Ao — A)^^^ 
is of class C^ in J7 X SI and satisfies the pointwise bounds 

^ WD^iAn - A)"^' h^^L^ \\D^{An - A)--'' U^^l^ 

«;c|Ao + Ar2j+5("+i"i+i^i). 

This estimate for the kernel of (Aq — A)^^-' readily yields the desired bounds for 
the Green's function. Indeed, the heat kernel H[t, x, y) (that is, the integral kernel 
of e*^*^) can be obtained from K{x^ y) as 



Hit,x,y)^ '-' ^^. J e-'^Kix,y)dX, 

where F is a contour enclosing the spectrum of Aq. Hence the C^ bound for K{x, y) 
readily implies that 

\DZD^yH{t,x,y)\ <;Ct-("+l"l+l^l)/2g(e-A^)t^ 

where e is an arbitrary positive constant. By Davies' perturbation method, this 
implies the Gaussian-type estimate 

\D'^DlH{t,x,y)\ s; c'i-("+l"l+l^l)/2e(^-^n)*e-^l^-^l'/* 

for i > 0. Finally, the estimate for the Green's function Gn{x,y) of the operator 
Ao — fc^ follows by elementary methods from the heat kernel's Gaussian bound and 
the identity 



/•OO 

Gn(x,2/)= / e-"'"H{t,x,y)dt. 
Jo 



'0 

a 

In order to effectively apply the previous theorem to the study of tentacled 
hypersurfaces, we will need the lower bound for the eigenvalues that we shall prove 
in the following lemma, which will ensure the C^ exponential decay at infinity of 
the Green's function of a tentacled domain: 

Lemma 5.6. Let fl be either a tentacled domain in M" or a Riemannian product 
of the form A x (T, oo) where T G [—00,00) and A is a bounded domain o/R"~^ 
with smooth boundary. Then An is strictly positive. 

Proof. It is a simple matter to show that Aax(t,oo) equals the lowest Dirichlet 
eigenvalue Aa of the bounded domain A C M"^^, which is positive. Hence let us now 
assume that fi is a tentacled domain. Since cannot be a Dirichlet eigenvalue of $7, 
the result will follow once we show that does not belong to the essential spectrum 
CToss(— An) of the Laplacian on 17 with Dirichlet boundary conditions. However, it is 
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well known that cross (—Aji) = acss{~^n\K) for any relatively compact set K C M" 
with smooth boundary. Therefore, taking K :— ilf] H^ n • • • n Hj , in the notation 
of Definition I5.1[ we find that 

JV 
CTossl-An) = CTcssl-AyN ^_5j — yj Crcss( — AsJ . 

i=l 

As the tentacle Si is isometric to the product A^ x (0, oo), it follows from our first 
observation that each Xs^ is positive, thus completing the proof of the lemma. D 

Before stating this section's main theorem, we need to establish some Green's 
function estimates for later use, the basic philosophy of which is to compare the 
Green's function Gq,{x, y) of a tentacled domain with that of a suitable domain with 
an Euclidean symmetry. To begin with, in the following lemma we will compare 
Gn{x,y) with the Green's function Gsi{x,y) of the tentacle Si when the points x 
and y belong to this tentacle: 

Lemma 5.7. Let Q be a tentacled domain. Then the pointwise estimate 

\D^D^Gnix,y) - D^DP^Gs,ix,y)\ ^ Ci/„^(x,y)e-^^['^-*(-^^-)+<^-*(^-'^-)l 

holds for all \a\,\/3\ ^ 2 whenever both points x and y lie in the tentacle Si. Here 
Ci and C2 are positive constants and we have set 

U^{x, y) := min { dist(x, A,)^"""!"! dist(y, A,)'-"-l^l , 

dist(a;, AO'-"-l"l dist(y, Kf-^'^P^ } . 

Proof. Let us take two distinct points x, y in the tentacle Si and apply Green's 
identity to Gsii'^x) and Gii{-,y) in this tentacle to derive the expression 



(5.4) Gn{x,y)-Gs,{x,y)= Gn{z,y)vi{z) -V ^GsAz,x)da{z) . 

Here da stands for the induced hypersurface measure on the 'tentacle cap' A;, Vi 
is the outer unit normal at A^ and we have used that G5. (-,2;) = on dSi and 

Gni-,y)=OondSAU 

Taking derivatives with respect to x and y in the identity (|5.4p and using Theo- 
rem [521 one readily finds that 



\D:D^Gn{x,y)-D:D^yGsAx,y)\^ f \D^Gn{z,y)\\V.D^GsAz,x)\da{ 

JAi 

^ ^3 J^^ |z-y|»+l/^h2|^_^|„+|.|-l '^^(^) 

C-f^QCile-iXn + k^y^^] dist(y,Ai)-|-C4[e-(Asi+fe^)'^^] dist(a;,AO 

^ dist(y, Ai)"+l/3|-2 dist(a;, A,)"+l"l-i 

for |a|, |/3| s; 2. As C2 := C4 min{(An + fc2)i/2 _ g^ (^Xs, + k'^f'^ - e} can be taken 
positive even if A; = by Lemma l5.6[ the above inequality proves the lemma with 

/„Mx, y) = dist(a;, A,)^-"-!"! dist(2/, A,)'-"-l^l . 

To show that this estimate also holds with 

fc.p{x, y) = dist(a;, A,)'-"-l"l dist(2/, Kf-^-\P\ , 
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thereby completing the proof of the lemma, it suffices to exchange x and y in 
Eq. (15. 4p by the symmetry of the Green's functions and repeat the argument. D 

In the following lemma we will prove the exponential decay of the Green's func- 
tion Gnix, y) when the points x and y lie in distinct tentacles: 

Lemma 5.8. Let x E Si and y G Sj be points lying in distinct tentacles of a 
tentacled domain fJ. Then the Green's Junction Gn{x,y) decays as 

\D^D^Gnix, y)| s^ Ci dist(a;, Ai)i-"-l"l dist(2/, Ai)2""-l'3le-^2[dist(a.,A.)+dist(y,A.)l 
for some positive constants Ci,C2 and all \a\, |/3| ^ 2 . 

Proof. Given x G Si and y € Sj, one can apply Green's identity to Gsi{-,x) and 
Gn{-,y) in Si to obtain the formula 



Gn{x,y) ^ / Gn{z,y)iyi{z) ■\'zGs,{z,x)da{z) . 
Ja, 

Taking derivatives with respect to x and y and using Theorem [531 £^i^ argument as 
in Lemma 15.71 vields 

\D^D^Gn(x,y)\^ f \D^Gn{z,y)\\V,D^GsA^,x)\da(z) 
Ja, 

r gC4[£~(An + fe2)i/2]|2_y| + c'4[e-(As.+fe^)^/^]|z-x| 

^ ^3 J^^ |z-y|»+l/^l-2|z-a;|"+l"l-i '^''^^^ 

C-i^eC4e-(\a+k^y^^ dist(y,Ai)]+C4[e-(As^+*;^)i/^] dist(i:,A.) 

^ dist(x, A,)"+l"l-i dist(y, A,)"+l'3|-2 ' 

The claim then follows by noticing that C2 := C4min{(Ao + fc^)^/^ — e, (A5. + 
/j2')i/2 _ g| ^^^ j-jg taken positive even for k — hy Lemma [5^ □ 

We shall next prove the main result of this section, namely, that any tentacled 
hypersurface can be transformed by a small diffeomorphism into a level set of a 
global solution to the equation (A — fc^)u = 0. The proof follows the strategy 
we outlined in Section [51 using the above estimates for the Green's functions to 
ensure that one can define a local solution of the equation that has a level set 
diffeomorphic to the tentacled hypersurface and satisfies the hypotheses of the C^ 
stability theorem: 

Theorem 5.9. Let L C M" be a (possibly disconnected) tentacled hypersurface of 
finite type and let k be a real constant. Then one can transform the hypersurface L 
by a diffeomorphism $ o/M", as close to the identity as we wish in the C^ norm, 
so that $(£) is a union of connected components of a level set m^^(co) of a solution 
of the equation (A — fc^)w = in R". 

Proof. We start by showing that, given any connected component Lq of the hyper- 
surface L, there exists a local solution of the equation having a level set diffeomor- 
phic to this component. To this end, let us denote by fl the tentacled domain whose 
boundary is Lq. We will keep the notation Si (with 1 ^ z ^ TV) for the tentacles of 
the domain H. (as in Definition 15.11) , which can be characterized as 

(5.5) S^■.= {y + tur■yeA^,t>0}. 
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Figure 3. Definition of tlie lialf-iine ji and the neighbor lioods V 
and U corresponding to a tentacie Si. 



Here the constant vector Vi is the outer unit normal at the tentacle cap A^. 

The construction of the desired local solution will make use of the Green's func- 
tion GQ{x,y) of the domain il. To ensure that the local solution satisfies the 
hypotheses of the C^ noncompact stability theorem, it is convenient to start by 
considering a straight half-line 7^ in each tentacle. That is, we fix a point yt in each 
tentacle cap A^ and define the half-line 7^ as 

(5.6) J, --^ {Vi + ti^i : t > 1} . 

The length measure on 7^ will be denoted by /i^ . A sketch of many of the geometric 
objects that appear in the proof of this theorem is given in Figure [S] 

The local solution will be constructed later on using the positive function 



Vt{x) := / Gn{x,y)d^ii{y). 

As a consequence of the estimates for the Green's function we proved in Theorem l5.5l 
and Lemma l5.6[ one can readily check that the function Vi satisfies 



Vi{x) ^ I Gn{x,yi+tvi)dt iS,C[ 

/•oo 

^C[ (dist(x,702+r2)' 



2;, -t;.,p-"e-^^l^-«' 



dt 



^-C'^(dmt(x.,iifWY''' 



dr 



— 00 
00 



/•OO 

< 2C[ / s2-"e-^^" ds 

J distix.- 



(5.7) 



€ 



t{x.,ti) 

'C^dist(a:,7,)^~"e'^2dist(a;,7,) jf „ ^ 4^ 



C^(l+ |l0gdist(x,7j)|)e-'^2dist(a;,7.)/2 if „ ^ 3 

for some positive constants C'. In the above inequalities, the variables t and r 
are related hy t := t — \yi — Yi{x)\, where Yi{x) is the unique point of 7^ such that 
|a; — yi(a::)| = dist(x, 7^), and s := (dist(a;, 7i)^-f t^)-'^/^. By the above estimates, Ui is 
well defined; indeed, it can be readily checked that it is of class C°°{^\^i)r\L\^^{Q) 
and satisfies the equation 

(fc^ - l\)vi = ^i 
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and the boundary condition Wili,, = 0. 

Our desired local solution will be the sum 

N 
i=l 

which is smooth in the closure of the domain Vt minus the union of all the half-lines 
7i and satisfies the equation 

N 

1=1 
with boundary condition v\lo — 0. Our next goal is to show that the function v 
satisfies the saturation, gradient and C^-boundedness conditions of the C^ non- 
compact stability theorem. That is, we claim that there exists a half-neighborhood 
V C ^ oi the component Lq and a positive constant 77 such that the function v 
satisfies the gradient condition 

(5.8) |Vw|^C>0 

in a set V that is saturated by v in the sense that any component of v^^{c) connected 
with V is contained in V for all < c < 77, and that moreover the second-order 
derivatives of the function v are bounded in V . 

In order to prove this claim, we introduce the auxiliary function 



(5.9) Vi{x) := I G§^{x,y)dfli{y) 

which is of class C°°{Si\^i) n L\ac{Si) by the same argument leading to the esti- 
mate (|5.7p . Here we are respectively denoting by 

S; := {y + ti^t-ye A,, t £ M} and 7, := {yi + ti^^:teR} 

the cylinder and straight line corresponding to the tentacle Si and to the half-line 
7i, while jli stands for the length measure on the line 7^. 

A simple symmetry argument shows that the value of the function Vi at an 
arbitrary point of the cylinder Si can be expressed in terms of the Green's function 
Ga- of the tentacle cap via 

(5.10) Vi{y + tiyi) = GA^{y,yi), 

where we are parametrizing the points in the cylinder as y -I- tvi as we did in 
Eq. (|5.5I) and we naturally identify the cap A^ with a bounded domain of IR"~^ 
using the property (iii) of Definition 15.11 The normal derivative of the function Vi 
at the boundary of the cylinder Si can similarly computed using the symmetry as 

(5.11) d^Vi{y + tvi)^duGA^{y,yi), 
where y is any point in the boundary of the cap A^ and t G M. 

By Hopf 's boundary point lemma \V6\ , it follows that the above normal derivative 
d,yG\. (•, yi) is strictly negative, so the boundedness of the cap Aj allows us to infer 
that 

|VGA.(-,y,)| ^Ci>0 
on a half- neighborhood Ui C Ai of the boundary dAi. Using again the fact that 
the cap Ai is bounded, it is standard that this set Ui can be safely assumed to be 
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saturated by the function GAi{-,yi), meaning that for any c at most one connected 
component of the level sets {y £ Ai : GAi{y,yi) — c} intersects Ui and that this 
component is actually contained in Ui. By the symmetry conditions (IS.lOp and 
(|5.1ip . this ensures that there exists a positive constant rji such that the gradient 
condition 

(5.12) |Vu,|>C2>0 

holds in the set Vi := {y + tvi : y S Ui, t £ M}. It should be noticed that all 
the derivatives of Vi are obviously bounded in Vi by symmetry. By the definition 
of the half- neighborhoods Ui, the set Vi contains a unique component of the level 
set v^^{c) for all values < c < 771. Therefore, the above discussion shows that 
the auxiliary function Vi satisfies the requirements of the C^ noncompact stability 
theorem. 

Motivated by this, our next step towards proving that the solution v also satisfies 
the conditions of the stability theorem is to control the difference between the local 
solution V and the auxiliary function Vi in the tentacle Si. To this end, let us take 
a point X € Si and estimate this difference as 



iD^vix) ~ D"u,{x)\ i^ / \D^Gn{x,y) ~ D^Gs,{x,y)\dn,{y) 
+ f \D^GsSx,y) - D:Gs^{x,y)\dfi,{y) 



(5.13) + Y. I \DZGn{x,y)\d^x,{y)+ [ \DZGs^{x,y)\dUy) ■ 

Since Si is a tentacled domain itself, one can now apply Lemmas 15.71 and 15.81 to 
obtain that for \a\ ^ 2 the first three integrals can be upper bounded by the 
exponential C3 Q-Cidist{x,Ki) ^ where Cj are positive constants. To estimate the last 
integral, let us denote by y^ ■= yt + Vi the endpoint of the half-line 7^ and apply 
Theorem 15.51 and Lemma 15.61 to derive that 

/■ Cc. p-C'ela;-!/! 

p:G,S^,y)\di,,{y)^ I d~^,[y)^C,e~^^\^-y^\ 

whenever the distance from the point x to the cap A^ is greater than 2 and \a\ ^ 2. 
Hence we obtain the pointwise C^ estimate 

(5.14) max \D°'v{x) - D°'v^{x)\ s$ Gg e''^^"^''^ , 

which holds in the tentacle Si provided that |a;| is large enough. 

Armed with these preliminary results, we can prove that the local solution v 
satisfies the hypotheses of the C^ noncompact stability theorem. As a first obser- 
vation, notice that, the domain n having a finite number of ends N, the gradient 
bound (|5.12l) and the G^ estimate (|5.14l) imply that there is a positive constant 772 
and a compact subset K of R" such that the gradient of the local solution satisfies 
\Vv\ ^ C" > in the set 

N 



W:^ [v-\{0,rj2))n\Jv)\K 
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and that the second-order derivatives of v are bounded in W. (We recaU each set Vi 
was chosen so that the auxihary function Vi and its derivatives satisfied appropriate 
bounds in it.) We can safely assume that, for any connected component Wq of W, 
there is a unique component of v~^{c)\K meeting Wq for all values c e (0, 7^2) and 
that this latter component of the level set does not intersect the boundary of Wq 
but at the compact set dK. 

The set W should be thought of as a conveniently chosen half-neighborhood of 
the hypersurface Lq minus a compact set. As the local solution v satisfies suitable 
gradient and saturation conditions in W by the above argument, now it essentially 
suffices to deal with v in the compact set K. Indeed, by Hopf's boundary point 
lemma [13] there are positive constants C2 and 773 and a half-neighborhood W C 51 
of the intersection Lq n K where the gradient of the local solution is bounded as 
\Vv\ ^ C2. Moreover, by compactness it is obvious that the set W can be chosen 
so that for all values c G (0,773) there is a unique component of v~-^{c) Ci K that 
meets the set W' , this component intersecting the boundary dW' only on dK. 

Putting together these results, Eq. ()5.8p . which ensures that the local solution 
V satisfies the gradient and saturation conditions of Theorem 13.11 now follows by 
taking the constant 77 := min{772, 773} and choosing an appropriate subset V C 
W U W . It also stems that the norm ||w||c2(v) is finite. 

Before we can profitably apply the stability theorem to the local solution w, 
there is one last technical point we must take care of. The equation {A — k'^)v = 
is only satisfied in the half- neighborhood V of the hypersurface Lq, not in its 
closure. Therefore, in order to apply the theorem we will first prove that there 
is a level set of u in V^ diffeomorphic to Lq via a C^-small diffeomorphism (e.g., 
via a diffeomorphism ^0 with \\'^q — id\\c^{R") < e/2) which only differs from the 
identity in a neighborhood of V. This is easily shown by taking an open set V 
containing the closure of V and a smooth extension v of the local solution v to 
the set {V U r2)\7 which is equal to v in i7\7 and negative in V'\i^ (but does not 
necessarily satisfy the equation (A — k'^)v = 0). One can then apply Corollary 13.31 
with {f,g,U,V,L) = {v,v - c,V' U (^^MJ^^i 7»), ^, ^o) and p = 1 to deduce the 
result, where c is a small enough constant in the interval (0, 7/). 

Applying the same reasoning to all the connected components of the hypersurface 
L, we derive that there exist a diffeomorphism ^I^ of R" with ||* — id||ci(Rn) < e/2 
and a function w, which satisfies the equation (A — k'^)w = in the closure of a 
neighborhood U of '^{L), such that: 

(i) The transformed hypersurface ^(L) is a level set w~^{cq) of the function, 

for some positive Cq. 
(ii) The neighborhood U is saturated, that is, if the intersection of w^^{c)nU is 
nonempty for some c G M, then w~^{c) does not intersect the boundary dU. 
(iii) The local solution satisfies the gradient condition |V7«[ ^ C" > in [/ and 

its second-order derivatives are bounded in this set. 
(iv) The complement of the set U in M" does not have any compact compo- 
nents. 

To complete the proof of the theorem, it suffices to approximate the local solution 
w in the set C/ by a global solution u of the equation (A — k'^)u = 0. By the 
condition (iv) above, one can invoke Theorem 14.51 to do so, ensuring that the C^ 
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norm |ju — u'||c'2(-[7-) is arbitrarily small. If this norm is chosen small enough, the 
construction of the local solution w and the saturated set U allows us to apply 
Theorem 13.11 to each connected component of the level set w~^{co) to obtain C^- 
small diffeoniorphisms that are only different from the identity in a prescribed 
neighborhood of the component and transform each component of w~^{co) into 
components of a level set of the global solution u. As there is no loss of generality in 
assuming that the supports of these diffeoniorphisms minus the identity are pairwise 
disjoint, we therefore obtain a diffeomorphism $ of K" with ||'I' — id||c'i(R") as small 
as we wish (say, smaller than e/2) transforming the level set w~^{co) into a union 
of components of a level set of u. The diffeomorphism $ := VJ/ o vl> then transforms 
the hypersurface L into a union of components of u^^(co) and is arbitrarily close 
to the identity in the sense that ||$ — idHcifRn-) < e. D 

Remark 5.10. The method of proof remains valid if we do not demand the ends 
of the tentacled domains to be 'straight' (i.e., isometric to Ai x (0, cxd)) but 'of 
Solomonic column type' (i.e., isometric to the intersection of A^ x (0,oo) with a 
domain invariant under a free isometric Z-action). 

6. NONCOMPACT LEVEL SETS: THE CASE OF INFINITE TOPOLOGICAL TYPE 

In this section we will conclude the proof of Theorem 1 1.1 1 bv considering the case 
of hypersurfaces that are not finitely generated. Although the basic philosophy 
of the proof is the same as in Theorem 15.91 in this case the hypersurfaces under 
consideration typically have an infinite number of ends (so, in particular, they 
are not diffeomorphic to an algebraic variety) that are not necessarily collared, 
and this introduces additional difficulties that require a separate treatment. The 
simplest example of a hypersurface of this kind is the torus of infinite genus in 
M.^ (cf. Figure [1]). This example shows that it is very convenient to embed the 
hypersurfaces so as to exploit discrete translational symmetries, so we will start by 
introducing some notation associated to these symmetry groups. 

Let us fix a positive integer i not greater than the space dimension n. We take 
a set of linearly independent vectors A :— {ai, . . . ,ae} C M" and denote by a* 
their dual vectors, which are the only elements in the linear span of the vectors A 
satisfying Ui ■ a* — Sij . For each i £ Z^ we will then denote by r/^ : R" — >■ R" the 
map 

T^{x) := X + hai H + tgag , 

which defines a free isometric Z^-action. We will also consider the fundamental cell 

Q-^ := {siai H 1- s^ae + b -.0 < Si <1, be M" and 6 • a^ = for 1 ^ i ^ £} 

associated to this action and the faces 

nf:^{xedQ-^:x-a*=0}, 
with li^ji^i. 

We will say a set U of R" is A-periodic if it is invariant under the above Z^ 
action, i.e., if t^{U) — U for all t &l/. If a set W is ^-periodic, U can be recovered 
from its intersection with the fundamental cell via the identity 

(6.1) U=\jT^{Ur\Q^). 
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For simplicity, we shall sometimes say that a set is i-periodic if it is ^-periodic for 
some set A with i independent vectors. 

Basically, the motivation of this section is to prove that there are global solutions 
to the equation (A — k^)u = having a level set diffeomorphic to infinite connected 
sums of any algebraic hypersurface. From the experience of Theorem 15.91 one can 
guess that it will be useful to exploit this diffeomorphism to embed the infinite-type 
hypersurface (in this case, the aforementioned infinite sum) so that both the collared 
ends of the underlying algebraic hypersurface and the way the different individual 
hypersurfaces are glued together are 'geometrically controlled' at infinity. We will 
do this through the following definition: 

Definition 6.1. An ^-periodic domain U of M" with smooth boundary is tentacled 
if its intersection with the fundamental cell Q is either relatively compact or equal 
to a tentacled domain of finite type J7 minus a compact subset K of M". A tentacled 
hypersurface of M" of possibly infinite type is a hypersurface with a finite number 
of connected components, each of which is the boundary of a (possibly periodic) 
tentacled domain. 

Remark 6.2. A tentacled hypersurface of infinite type does not need to be periodic, 
even if all its components are. Moreover, the periodic components can have distinct 
symmetry groups. 

It should be noted that if the intersection of the periodic tentacled domain U 
with the fundamental cell Q is unbounded, obviously the rank £ of the symmetry 
group is at most n — 1. 

The class of tentacled hypersurfaces of infinite type modulo diffeomorphism in- 
cludes infinite connected sums of a large class of hypersurfaces, as we will see in 
the following examples: 

Example 6.3. If L is a (possibly compact) algebraic hypersurface of M", there is 
an .^-periodic tentacled hypersurface that is diffeomorphic to a connected sum of 
infinitely many copies of L. Here the rank £ can take any value between 1 and n — 1 
(resp. n) if L is noncompact (resp. compact). In particular, and getting back to 
Question [TJ the torus of infinite genus and the infinite jungle gym are examples of 
1-periodic and 3-periodic tentacled surfaces of M^, respectively. 

Example 6.4. Given any integer £ between 1 and n—1 and a tentacled hypersurface 
L in R", there is an ^-periodic hypersurface that is diffeomorphic to a connected 
sum of infinitely many copies of L. This readily follows from the following ele- 
mentary proposition, which is a trivial consequence of the fact that any tentacled 
hypersurface is collared: 

Proposition 6.5. Given a tentacled hypersurface L C K" and a set A C K" of £ 
independent vectors (£ !^ n — 1), one can transform it by a diffeomorphism ^ o/M" 
so that ^(£) is tentacled and contained in the fundamental cell Q . 

As in the previous section (cf. Lemma 15.61) , firstly we need to prove that the 
spectrum of the Dirichlet Laplacian in a periodic tentacled domain is bounded away 
from in order to obtain exponential decay of the Green's function. This is done 
in the following lemma, which exploits both the asymptotic Euclidean symmetries 
of tentacled domains of finite type and the invariance under the isometric Z^ action 
of periodic tentacled domains: 
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Lemma 6.6. LetU he an A-periodic tentacled domain in M". Then the bottom of 
the spectrum of the Laplacian Xjj in this domain is positive. 

Proof. Let us begin by observing that the Dirichlet spectrum in U can be written 
as 

(6.2) ct{Au)^ \J a{Ae), 

where cr{Ag) denotes the L^ spectrum of the Laplacian in WflQ'^ with the boundary 
conditions 

(6.3) w = oiidUnQ^, 

D°'w{x + Oj) = e''^W^w{x) for al\x eUD Uf, \a\ sC 1, 1 ^ j s^ ^. 

Here T'^ := (R/2ttZY is the £-torus and we write = {0u . . . ,9e). As for Eq. (jg:^ . 
the inclusion of a{Au) in the union of the spectra of A^ follows from a standard 
modification of Floquet theory [5S] , while the fact that both sets actually coincide 
follows e.g. from Adachi's results on unitary actions of amenable groups [2]. 

From the boundary condition (|6.3|) and the fact that the boundary dU intersects 
the fundamental cell Q-^, it follows that cannot be an eigenvalue of Ag for any 
6 E T^. Besides, cannot belong to the essential spectrum of Ag either, since 
the spectrum ^(Ajypqa) is bounded away from and (JcsslAe) coincides with the 
Dirichlet essential spectrum aessiAnf-^Q^) by the boundedness of the intersection 
U n n^. The former assertion is clear when VI fl Q-^ is bounded and stems from 
Lemma [5.61 when U n Q"^ is a tentacled hypersurface minus a compact set. The 
statement now follows from the decomposition (16. 2p and the compactness of T^ . D 



We shall next prove the main result of this section where we adapt the method 
of proof of Theorem l5.9l to construct solutions of the equation {A — k'^)u = with a 
level set diffeomorphic to any tentacled hypersurface of infinite type. To avoid un- 
necessary repetitions, we will not present in full detail some steps in the argument, 
referring instead to the appropriate parts of the demonstration of Theorem 15.91 
Together with Theorem 15. 9[ this completes the proof of Theorem 11.11 

Theorem 6.7. Let k be a real constant. Given a (possibly disconnected and of 
infinite type) tentacled hypersurface L C M", we can transform it by a diffeomor- 
phism $ o/R", arbitrarily close to the identity in the C^ norm, so that ^{L) is a 
union of connected components of a level set u^^(cq) of a solution to the equation 
{A-k^)u = inM". 

Proof. As in the proof of Thcorem l5.91 our goal is to construct a local solution of the 
equation defined in a half-neighborhood of each component Lq of the hypersurface 
L and satisfying the saturation, gradient and C^-boundedness conditions of the C^ 
noncompact stability theorem. If the component Lq is a tentacled hypersurface 
of finite type, this local solution was constructed in the proof of Theorem 15. 9i so 
we will assume that Lq is y^-periodic for a set of £ linearly independent vectors A. 
We will denote by U the ^-periodic tentacled domain enclosed by Lq; the main 
geometric objects considered in this proof are presented in Figure |31 

By the definition of an ^-periodic tentacled domain, the intersection 14 n Q-^ of 
this domain with the fundamental cell is either relatively compact or equal to ^\K, 
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7l 



72 
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Figure 4. The half-lines 7i, the fundamental cell Q-^ and the 
neighborhood V that appear in the proof of Theorem 16.71 



with O a tentacled domain with N ends and K a compact set. Let us first suppose 
that U n Q-^ = Q\K. In this case, we can safely assume that the tentacles Si of 
the domain f2 do not intersect the compact set K, and define the half-lines "fi C Si 
as we did in Eq. (|5.6|) . Let us set 



Vt{x) := / Gu{x,y)dm{y), 

with /ii being the length measure on the half-line 7^ . By Lemma 16.61 and the 
exponential bound for the Green's function proved in Theorem l5.5[ one immediately 
obtains the estimate 



Ci dist(a;, 7,)3-"e^'^2 dist(:r,7.) if n ^ 4, 

Ci(l + |logdist(a;,7,)|)e-^2'^''''("^'T-)/2 jf ^j ^ 3 



Vi{x) < 



for the function Vi after arguing as in Eq. (|5.7|) . 

The symmetry of the domain U implies that its Green's function is invariant 
under this symmetry group: Gi{{x,y) = Gu{T^{x),T^{y)). Therefore, the above 
estimate implies that the function 

N 

is well defined, for the above sum converges uniformly on compact subsets of U\"f, 
with 

7 := U U ^t-'M 

being the union of all suitably translated copies of the half-lines 7^ . The function 
V, which is of class C°°{U\j) fl L\^^{U), therefore satisfies the equation 



N 



ie-A)v = j2Y.i-t% 



fj^i 



i=l t€ 



in U and the boundary condition v\lo — 0; indeed, we shall next show that v is the 
local solution we wanted to construct. In the above equation, (t^)* denotes the 
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push-forward. From its definition it is clear that v inherits the symmetries of the 
domain, meaning that 

(6.4) v{t^{x))^v{x). 

The next step of the proof consists in showing that the local solution v satisfies 
the hypotheses of the C^ noncompact stability theorem, that is, there are a half- 
neighborhood V C U oi the component Lq with ||u||c2(V') < cxd and some positive 
constant r/ such that miy |Vw| > and any component of w~^(c) connected with V 
is contained in V for all < c < r/. By the symmetry of v (cf. Eq. (j6.4[) '). it suffices 
to prove the claim in Z// fl Q^. Since U Ci Q-^ equals il\K, a close look at the proof of 
Theorem 15.91 shows that the proof we gave for tentacled hypersurfaces carries over 
almost verbatim to the present situation. Indeed, the difference between D°'v{x) 
and the function D"vi{x) introduced in the relation (|5.9|) can be computed as in 
Eq. (|5.13p . the only changes being that one must use the Green's function Gu{x, y) 
instead of Gn(x, y) and the sum over the set of integers {j : 1 ^ j ^ N and j ^ i\ 
is to be replaced by a sum over 

{(j, t):li^js^N, t e Z^ and j / i if i = 0} . 

As Lemmas 15. 71 and 15. 81 and their proofs obviously remain valid with the domain Q, 
replaced by U^ on account of Lemma l6.6l the rest of the argument remains unchanged 
and the above claim follows. This allows us to apply Corollarv l3.3l as in Theorem l5.9l 
to show that the level set v~^{cq) is difFeomorphic to the component Lq for any 
sufficiently small positive constant cq. 

Let us now treat the case when the intersection U n Q^ is relatively compact. 
Then we take a point y^ in this intersection and let 

vo{x) := Gu{x,yo) 

be the Green's function of the domain lA with pole at yo- By Lemma 16.61 and 
Theorem 15.51 one has the exponential bound 

vo{x)^Ci\x-yQ\'-'\-^-\^-y°\, 

which allows us to show that the function 

v{x) := Y,vo{t^{x)) 



is a well defined solution of the equation 

(fc^ - A). = ^ 5,.(,,) 

with the symmetry property (|6.4p . Therefore to show that the function v is the 
local solution we are looking for it is enough to prove the gradient bound |Vf| ^ 
C > and the boundedness of the second-order derivatives of i; in a saturated 
neighborhood of the set Lq n Q^^^ which is straightforward by Hopf's boundary 
point lemma and the compactness of Lq n Q-^. The rest of the argument goes as in 
the previous case, where U n Q^^ = ^\K. 

Hence, applying the same reasoning to all the connected components of the 
hypersurface L we can use the same argument as in Theorem 15.91 to infer that 
there is a diffeomorphism ^ of M" with ||^ — id||c'i(Rn) < e/2, a neighborhood U 
of ^{L) and a local solution w of the equation (A — fc^)z« = in the closure of U 
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which satisfy the properties (i)-(iv) in the proof of Theorem lS.QI The theorem then 
follows from the Approximation and Noncompact Stability Theorems 14.51 and 13.1 
as in the proof of Theorem 15.91 D 



Remark 6.8. It is clear that the result and the method of proof remain valid also 
if L is allowed to have an infinite number of connected components provided we 
impose appropriate 'uniform' assumptions. For example, it suffices to impose that 
the components of L do not accumulate and that each component be isometric to an 
element of a fixed, finite collection {Li, . . . ,Lr} of tentacled hypersurfaces (possibly 
disconnected and of infinite type) . Likewise, the method of proof also works when 
the domain U is not necessarily periodic but there is a periodic tentacled domain 
U and a compact subset K of R" such that U\K = IA\K . 

7. Joint level sets 

In this section we shall study joint level sets of solutions to elliptic equations, 
which amounts to considering transverse intersections of level sets. In doing so, 
we will find it convenient to consider separately the cases where the components 
of the joint level set are all compact or not. This difference in treatment is due to 
the fact that when all the components are compact one can exploit a fairly general 
better-than-uniform approximation result to deal with a wide class of equations, 
while when there are some noncompact components the situation is considerably 
more subtle and our approach relies on our results for single level sets. (It is worth 
analyzing why the treatment of a single compact level set, for equations that admit 
these kind of level sets, is considerably less involved than the case of noncompact 
level sets; for this, we refer to Appendix \^ where this problem is considered in 
detail.) 

7.1. The compact case. Here we shall consider the compact joint level sets 
u^ (ci) n • • • n u'^{cm) of solutions to equations 

Tr-Ur — 

in M". Here each Tr is a linear elliptic differential operator of second order with 
real analytic coefficients and the number m of functions is at least 2. As we shall 
see, the reason why we can consider more general equations in this section is that 
we can prove a better-than-uniform approximation theorem for locally finite unions 
of disjoint compact subsets of R". This theorem provides fine control at infinity, 
which can be combined with the Cauchy-Kowalewski theorem in each component to 
derive the desired realization results. This approach does not work when some of the 
components are noncompact because the approximation result given in Theorem l4.5l 
is not fine enough to deal with the domains of definition of the local Cauchy- 
Kowalewski solutions, which can be very narrow at infinity. 

For the sake of completeness, we start with the following standard lemma, which 
allows us to approximate a smooth submanifold by an analytic one (we will later 
on apply this result in the case of codimension 1 in the preparations to apply the 
Cauchy-Kowalewski theorem) : 

Lemma 7.1. Let L he a locally finite union of pairwise disjoint compact codimension- 
m suhm,anifolds of M" with trivial norm,al bundle. Then one can transform L by a 
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dijfeomorphism ^ o/R", arbitrarily close to the identity in any C^ norm, so that 
^(L) is an analytic submanifold. 

Proof. Let Li, denote the connected components of the submanifold L, with b rang- 
ing over an at most countable set B. Since Li, has trivial normal bundle, there is 
a tubular neighborhood Wb of Lb and a C°° trivializing map Qb ■ Wb -^ K™ with 
6^^(0) = Lb- We can safely assume that the closures of the sets Wb are pairwise 
disjoint. For any 5b > 0, Whitney's approximation theorem [501 Theorem 1.6.5] 
enables us to take a real analytic submersion 9(, : Wb -^ R™ with ||6b — QbWcpiWb)- 
The hypotheses of Theorem 13.11 being automatically satisfied by the compactness 
of the component Lb and the fact that the map Qb is a trivialization, we can now 
apply Theorem l3.1l and Remark [3^ with (/, g, U, V, L) equal to (9b, Qb, Wb, Wb, Lb) 
to derive the existence of diffcomorphisms ^I^ ^ of M" such that: 

(i) Q^\Q) = ^b{Lb). 
(ii) ll^b — id||cp(R") < e and the support of '^b — id is contained in Wb- 

Hence these diffeomorphisms naturally define the desired diffeomorphism \1' by 






D 



Instead of Theorem 14.51 in this section we will use the following better-than- 
uniform approximation result, which is modeled upon a theorem of Bagby and 
Gauthier [4]. This result, which is valid for general analytic elliptic operators 
but requires an essential compactness assumption, is proved through an iterative 
argument using the Lax-Malgrange theorem (which accounts for the conditions we 
impose on the set S) and a suitably chosen exhaustion of M": 

Lemma 7.2. Suppose that S C M" is a locally finite union of compact sets of 
nonempty interior whose complements do not have any relatively compact compo- 
nents. If w is a local solution of the equation T^-w = in the set S, then it can 
be approximated in S by a global solution of this equation in the C^ better-than- 
uniform sense. That is, for any p and any positive continuous function e{x) there 
is a function v satisfying the equation T^v = in R" such that, pointwise in S , 

max \D°'v{x) ~ D°'w{x)\ < e{x) . 

\a\^p 

Proof. To control the C^ norm of the difference v — w as required by the function 
e(x), we will introduce some positive constants e^ associated to this function and 
to an exhaustion of R" that we shall define next. For this, let us denote by Sb the 
connected components of the set S, where b ranges over an at most countable set 
B. As S" is the locally finite union of Sb, we can take an exhaustion =: Kq C 
Ki C K2 C ■ • • of R" by compact sets such that: 

o 

(i) The union of the interiors Kj of the sets Kj is the whole R" . 
(ii) For each j, the complements of Kj and of 5 U Kj are connected. 
(iii) If some Sb intersects Kj, then Sb is contained in the interior of Kj+i. 
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Now we can take any positive numbers Cj such that 



(7.1) 



^ oo 

e., < — min e(x) and > £*■ < e 



k=J + l 



for all j ^ 1. We also set eq '■— 0. 

We now proceed by induction. We make the induction hypothesis that there are 
functions Vj satisfying T^-Vj = in M" and such that, for any s ^ 1, 



(7.2a) 

(7.2b) 
(7.2c) 



S 



< e,s 



CP{Sn(K, + i\K,)) 



<es + 2es 



CP{Sn{K^\K,-i)) 



McHK.-i) <^s + es-i. 



Let us begin by proving the induction hypotheses for s = 1. As the complement 
of SnK2 does not have any relatively compact components, the Lax-Malgrange the- 
orem [301 Theorem 3.10.7] yields a function vi which satisfies the equation T^ui =0 
in M" and such that 

lh-^i|lcp(SnK.)<^i- 

By the definition of the set Kq and of the constant eo , it is therefore evident that 
the induction hypotheses hold in this case. 

Let us now assume that the induction hypotheses hold for all 1 ^ s ^ fc and 
use this assumption to prove them for s = k + 1. For this purpose, let us define a 
function Wk on 5 U Kk by setting WklKk '■— and 



w - J2 ""j if Sb intersects Kk+2\ Kk+i, 

Wk\St ■= \ = 1 





if Sb does not intersect Kk+2\ Kk+i ■ 



The definition of the exhaustion and the hypothesis (|7.2ap guarantee that TrWk = 
and 



(7.3) 



Ffe|lcp(KfcU(SnKfc+i)) ^ 



fc 



< Cfe- 



Cp(sn(Kfc+i\Kfc)) 



Since the complement of the set Kk+2 H (5 U Kk) does not have any bounded 
components by the conditions we imposed on the exhaustion, a further application 
of the Lax-Malgrange theorem allows us to take a solution Vk+i of TrVk+i = in 
M" such that 



(7.4) 



\wk - «'c+i||cp(K,+,n(suK,)) < ^^=+1 ■ 



Eq. (|7.4p and the definition of Wk ensure that the hypothesis (j7.2ap also holds 
for s = fc + 1. Moreover, Eqs. ()7.2ap . (|7.3p and \7 A\ imply the following pointwise 
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D-iw-Y^. 



CP bound in the set 5 n {Kk+i\Kk), valid for all \a\ < p: 

3 = 1 ' ^ J = l ^ 

< eA; + \D°-{wk -Wfc+i)| + |£'"«^fe| < efc+i+2efc, 
This proves the second induction hypothesis (j7.2bl) for s = fc + 1. Moreover, 

lhfc+i|lcp(Kfc) ^ lhfc-^fc+i|lcp(Kfc) + lkfc|lcp(i<-,.) <efc+i+efc 

as a consequence of Eqs. (I7.3P and (j7.4p . so the remaining induction hypothe- 
sis (|7.2cl) also holds for s = fc + 1. The induction argument is then complete. 

Let us now define the global solution v as 

the sum converging C-uniformly by the definition of constants tj and the induction 
hypothesis (|7.2cp . Since each Vj satisfies TrVj = 0, this C convergence ensures that 
the function v satisfies the equation T^w — too. Besides, from the conditions (|7.ip 
we imposed on the constants tj and the induction hypotheses (|7.2p it follows that 
in each set S n {Kk+i\Kk) one has the pointwise C estimate 



\D"'{w-v)\ < 



fc+i 



^"h-E"^ 



j=i 



li^^i'M 



V 4 — h^-i ' 



j=k+3 



< (efc+i + 2efe) + (efc+2 + ffc+i) + E ('^■J' + ^^-i) 

i=fc+3 

< 2efe + 4efc+i < min e(x) 

xG-fffc+i 

for any fc and all \a\ ^ p, as we wanted to show. 

We now have all the ingredients needed for the proof of Theorem 11.31 



D 



Proof of Theorem \1.3\ Our goal is to obtain the local solution by means of a Cauchy 
problem. For this, it is convenient to observe that Lemma 17.11 ensures that, by 
perturbing the submanifold a little if necessary, there is no loss of generality in 
assuming that L is a real analytic submanifold of R" . 

Let us denote by Lf, the connected components of L, with b taking values in 
an at most countable set B. We start by realizing Lf, as the intersection of m 
hypersurfaces Eb^- As each component Lf, also has trivial normal bundle, we can 
take an analytic trivialization 8;, : Wb -^ M™, where Wb is a tubular neighborhood 
of Lb and 6^^(0) = Lf,. We denote the components of 6b by {6bi, ■ ■ ■ ,Obm) and 
consider the analytic hypersurfaces Y,br '■— ^br^(O) C Wb, with 1 ^ r ^ m. By the 
definition of Ob, it is apparent that these hypersurfaces intersect transversally at 
Lb = Sbi n • • • n Tibm- 

Now that we have expressed the component Lb as the intersection of m real 
analytic hypersurfaces Sbr, we can consider the following Cauchy problems: 



TrVbr = , 



Vbr\s, 



duVb 



rISb 



= 1. 



40 ALBERTO ENCISO AND DANIEL PERALTA-SALAS 

Here du denotes a normal derivative at the hypersurface Yibr- The differential 
operator Tr being analytic and elliptic, the Cauchy-Kowalewski theorem grants the 
existence of a solution Vr to this Cauchy problem in the closure of a neighborhood 
Uhr of the hypersurface Eb^- As the hypersurfaces Sbr intersect transversally at Lb 
and the gradient of the solution vi^r coincides with the unit normal of Yii,r on this 
hypersurface, it stems that 

Ch '■— inf min wi \/vbi{x) + • • • + uj„i Vvbrnix)] 

is positive. By the compactness of Lb, the sets Ubr can be chosen small enough so 
that 

(7.5) inf min Li Vubi(a;) H h w™ Vubm (a;) I > -— , 

x&Jb |cj| = l Z 

where the set Ub '■= Ubi n • • • fl Ubrn can be assumed to be saturated by the functions 
Wbi , ■ . ■ , Vbm without loss of generality. We can also suppose that the sets Ub are 
pairwise disjoint. 

Let us set U := Ubes ^b and define the local solutions Vr of the equation TrVr — 
on the set U as Vr\ui, '■— "^br- By Lemma [7.21 this local solution can be approximated 
in the C^ better-than-uniform sense by global solutions Ur of the equations TrUr = 
0. That is, for any given positive continuous function e{x) in U (to be specified 
later) we can assume that 

max \D°'Ur{x) — D°'Vr{x)\ < e{x) 



in the closure of the set U. In view of Eq. (|7.5|) and the compactness of Ub, 
Theorem 13.11 and Remark 13.21 guarantee the existence of a constant eb > such 
that if e{x) < Cb in Ub there is a diffeomorphism <i>b of M", arbitrarily close to the 
identity in the C^ norm and equal to the identity outside Ub, which transforms the 
submanifold Lb into a connected component of the joint level set Ui (ci) fl • ■ • n 
M~^(cm). As the supports of each map <i>b — id are pairwise disjoint and we can 
choose £{x) so that it is smaller than eb in each set Ub, the theorem then follows 
by letting the diffeomorphism $ be equal to <l>b in each set Ub and equal to the 
identity in the complement of [/ = IJbes ^b- D 

Example 7.3. It is known that any exotic ?Ti-sphere smoothly embeds in M'^™ with 
trivial normal bundle |14| and that any locally finite link in R"^ has trivial normal 
bundle too, so Theorem 11.31 obviously furnishes a positive answer to Questions [2] 
and [31 In particular, there are two harmonic functions ui,U2 in K"^ such that 
u^^{0) n U2^{0) contains a knot in each isotopy class. 

Remark 7.4. When L is a finite union of compact codimension-m components and 
Tr = A for all r, one can actually prove that ui, . . . ,Um can be chosen to be 
harmonic polynomials by proceeding as above and using an approximation theorem 
of Paramonov [21, instead of the Lax-Malgrange theorem. 

7.2. The noncompact case. In this subsection we present a realization theorem 
for noncompact joint level sets u^ (ci)n- • ■nu^{cm) of solutions Ur of the equations 
{A — kf.)ur = 0, for any real constants fc^. It is assumed that the number of functions 
m is at least 2. As in the case of a single level set, it will be crucial to control the 
geometry of the submanifolds at infinity using asymptotic translation symmetries. 
In order to implement this idea, we will resort to the codimension-TTi analog of a 
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tentacled hypersurface, which can be defined in terms of tentacled hypersurfaces in 
a very simple manner: 

Definition 7.5. A (periodic) codimension-m tentacled submanifold is the trans- 
verse intersection of m tentacled hypersurfaces (possibly disconnected and of infi- 
nite type). 

Before going on, let us briefiy discuss why this is the suitable codimension- 
m analog of a tentacled hypersurface. From Definition I5.1[ it is clear that the 
properties one would require for a submanifold L ( for concreteness, connected and 
of finite type) to be a codimension-m analog of a tentacled hypersurface are the 
following: 

(i) L must have trivial normal bundle, as this is a topological obstruction to 
being a transverse intersection of level sets. 

(ii) There must be J embedded images Ilj of K"~^ in M", which divide R" 
into two domains HJ^ and H~ satisfying certain properties analogous to 
the conditions (i)-(iv) in Definition 15.11 In particular, each component of 
the intersection of L with Hi U • • • UII,/ is compact, has dimension n — m~l 
and is contained in an affine {n — rn)-plane. If its components are called 
Ai , the submanifold L has a finite number of ends and are all isometric to 
some product A^ x (0,oo). 

That the codimension-?Ti submanifold L is then the transverse intersection of tenta- 
cled hypersurfaces L^, . . . , i™ readily follows from the fact that, as a consequence 
of these properties, L has a tubular neighborhood whose boundary is a tentacled 
hypersurface. A very similar argument is also valid in the periodic, disconnected 



Example 7.6. In particular, from Example 15.41 it follows that the transverse in- 
tersection of m algebraic open hypersurfaces is diffeomorphic to a codimension-m 
tentacled submanifold. For instance, the torus of genus g with N ends can be 
realized as a 2-dimensional tentacled submanifold in any R". 

We shall next present the proof of Theorem ll.4[ which asserts that if L is a tenta- 
cled submanifold of codimension ni and L is a finite union of compact submanifolds 
of the same codimension with trivial normal bundle, then there are m solutions u^. 
to the equations (A — fc^)ur = in R" having a joint level set diffeomorphic to iUi 
(provided, of course, that L and L are disjoint). It is apparent that the 'compact' 
part of the argument can be easily tackled arguing as in the proof of Theorem II. 3[ 
so most of the proof will be devoted to the case of tentacled submanifolds: 



Proof of Theorem \1.4\ Let Lq be a connected component of the codimension- ?7i 
submanifold LU L. We want to prove that there are m local solutions u,- of the 
equations (A — k'^)vr = 0, defined in a suitable domain whose closure contains Lqj 
which have a joint level set diffeomorphic to Lq and satisfy the hypotheses of the C^ 
Stability Theorem 13.11 When Lq is compact, this immediately follows from Theo- 
rem [131 so we will henceforth assume that Lq is noncompact. By hypothesis, Lq 
is then given by the transverse intersection of m (connected but possibly periodic) 
tentacled hypersurfaces L^, . . . , L™. 

Let us concentrate for the moment in the tentacled hypersurface L^. As we 
showed in the proof of Theorem 15.91 (or 16. 7| if L'' is periodic) , there is a local 
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solution Vr , defined in a half- neighborhood Vr of L"^ and satisfying the conditions 
(i)-(iv) in the proof of Theorem 15.91 namely: 

(i) There is some positive and arbitrarily small constant whose corresponding 
level set v~^(cr) is diffeomorphic to the hypersurface L"^. This difFeomor- 
phism can be chosen supported in a small neighborhood of Vr and close to 
the identity in the C^ norm, 
(ii) The neighborhood Vr is saturated by the function Vr- 

(iii) The gradient of the function Vr satisfies the lower bound \Vvr\ ^ C > 
in the set Vr and on Lq can be written as Vvr = — iVu^l ^^n where z/^ is 
the outer normal at the hypersurface. Besides, the C^ norm Hwrllc^cyr) ^^ 
finite. 

(iv) The complement of the set Vr does not have any bounded components. 

We claim that the joint level set w^^(ci) fl • • • fl v:^j^{cm) is diffeomorphic to the 
component Lq via a C^-small diffeomorphism. To prove this, we start by noticing 
that the set V := V^ifl- • -nVm is saturated by the map (ui, . . . , Vm), which moreover 
has bounded second-order derivatives in V by the condition (iii) above. To show 
that the gradient condition of the stability theorem holds, we aim to prove that for 
any unit vector w — (wi, . . . ,uJm) the function 

I 1 2 

Fuj{x) := \iOiVvi{x) -\ 1- WmVwm(a;)| 

satisfies the bound 

Fu,{x) ^ Ci 
in the set V for some positive constant Ci independent of w. As the hypersurfaces 
L^, . . . , L™ intersect transversally at the component Lq and their geometry is con- 
trolled at infinity because they are tentacled, an easy geometric argument shows 
that the quantity 

I 1 2 

C2 '■= inf min luji vi [x) + • • ■ + uJm I'mix) 

xELq \ui\=1 

is positive. As Vvr = —\SIvr\vr on Lq with |Vwr(a;)| 5^ C by the condition (iii) 
above, this immediately implies that, for any x ^ Lq, 

F^{x) = |a;i|Vwi(a;)| vi{x) H (- Wm|Vw„i(x)| Vm{x)\ ^ C^Ci 

is bounded from below by some positive constant independent of w. As the G^ 
norm of the functions Vr is bounded and their gradients satisfy \Vvr\ ^CmV, it 
is straightforward to check that, by taking the saturated set V smaller if necessary, 
we can assume that inf|^|^i mixev F^ix) is also positive. 

As a consequence of the above digression, we can apply Corollary 13.31 with 
(/,5,i) — ((wr)J^i,(wr — Cr)^^i,LQ) (thc domain V in Corollary 13.31 coincides 
with the current set V, while U is a. suitably chosen neighborhood of V). This 
way we infer that, for any sufficiently small positive constants Cr, the joint level 
set Vi (ci) n • ■ ■ n v^^{cm) is diffeomorphic to the component Lq via a C^-small 
diffeomorphism supported in a neighborhood of V, as we claimed. 

Applying the same argument to each connected component of the submanifold 
L we obtain local solutions Wr of the equations (A — k'^)wr = that have a joint 
level set diffeomorphic to L via a C^-small diffeomorphism of R" and satisfy the 
hypotheses of the C^ stability theorem. The number of connected components of 
L being finite, this allows us to apply the Approximation Theorem 14.51 and the 
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Stability Theorem 13. II to complete the proof of the theorem (the details follow just 
as in the proof of Theorem 15.91 and are thus omitted) . D 

Remark 7.7. Obviously, the analogs of Remarks 15. 101 and 16.81 also apply to the case 
of joint level sets. 

Appendix A. The case of a compact level set 

For the sake of completeness, in this Appendix we shall consider the construction 
of solutions to the equation (A — q)u = in M" having a level set diffcomorphic to 
a union of compact hypersurfaces. 

A compact hypersurface separates R" into an inner (bounded) domain and an 
outer (unbounded) domain. In the following theorem we shall assume that the inner 
domains associated to each component of the hypersurface are pairwise disjoint. 
The reason for this is that a key element of the proof is the better-than-uniform 
approximation result established in Lemma l7.2[ which makes use of this hypothesis. 
Upon comparison with the proof of Theorem 11.11 this illustrates how the proof 
drastically simplifies in the case of compact components, as many of the subtleties 
that appeared in the former case are now absent. 

Theorem A.l. Let L be a locally finite union of compact hypersurfaces o/M" and 
let q be a nonnegative, real analytic function on M" that is not identically zero. 
Suppose that the inner domains corresponding to the connected components of L 
are pairwise disjoint and fix a positive constant c. Then one can transform the 
hypersurface L by a diffeomorphism $ ofM." , arbitrarily close to the identity in any 
CP norm, so that ^(L) is a union of connected components of the level set u^^(c) 
of a function that satisfies the equation (A — q)u — in M". 

Proof. Let us denote the connected components of L by Li,, the index b ranging 
over an at most countable set i?, and let flf, stand for the inner domain of the 
component Li,. We saw in Lemma F7. II that, perturbing L a little if necessary via 
a C^-small diffeomorphism, there is no loss of generality in assuming that all the 
components Lf, are real analytic. 

A natural way of constructing a local solution of the equation having a level 
set diffcomorphic to Lf, is via a boundary value problem in the domain £7^. It is 
standard that the boundary value problem 

(A-(j)wb=0 in fib, Wb|Li,=c 

has a unique solution, which is given by 



Vb{x):=c-cl Gn,{x,y)q{y)dy. 

Here Gn^{x,y) is the Green's function of the operator A — g in the domain fib, 
defined as in Proposition l4.1l The function vj, satisfies the equation (A — q)wb = in 
the closure of J7f, by the fact that the function q and the domain fit are analytic 118] , 
is smaller than c in Vlb by the maximum principle, and its gradient Vuf, is nonzero 
on Lb by Hopf's boundary point lemma |13j . Hence we can now define a function 
V on the closed set S := Uie/ ^fc that satisfies (A — q)v = by setting v\ji := Vb. 
Obviously v^^{c) = L. 
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Figure 5. A solution of the equation Au = u in R'^ whose zero 
set contains a compact surface of each genus. 



We shall next prove that one can use Lemma 17.21 on better-than-uniform ap- 
proximation to find a global solution u of the equation such that L is diffeomorphic 
to a union of connected components of u~^{c). This makes use of Remark 13.21 on 
the topological stability theorem, from which it stems that, given any component 
Lb and any positive constant ep, there is a positive constant Sf, such that for any 
function u with ||m — w||cj'(Oi,) < ^b one can find a diffeomorphism $b of R" with 
W^b — id||cp(R") < £o mapping Lf, onto a connected component of u~^(c). It should 
be noticed that the saturation and gradient conditions of the stability theorem are 
obviously satisfied as a consequence of the gradient estimate Vwb 7^ on i;, and 
the compactness of each component. One can also assume that $b — id is supported 
in a neighborhood Ub of fib, with the sets Ub pairwise disjoint (cf. Theorem 13. 1|) . 

To complete the proof of the theorem, let us take a positive continuous function 
e{x) which is smaller than Sb in each set f2f,. By Lemma 17.21 there is a function u 
satisfying the equation {A — q)u = in R" and the following C^ better-than-uniform 
bound in the set S: 

max iD^uix) ~ D°'v{x)\ < €{x) . 

By construction, $(i) is then a union of connected components of u^^(c), where 
the diffeomorphism $ is defined as ^{x) := $b(x) if x belongs to some set Ub and 
^{x) := X otherwise. D 

Remark A. 2. It is clear that the same proof works for a wider class of positive 
second-order operators, but we will not purse this issue here. The case of the 
equation (A + k'^)u = 0, considered in Question [5] in the context of joint level 
sets, can also be easily dealt with. To construct the local solution having a zero 
set diffeomorphic to L, it suffices to consider the first Dirichlet eigenfunction Vb of 
each domain fib, which satisfies (A -|- Xb)vb = for some constant Xb > 0, and to 
apply suitable dilations and rigid motions to the functions Vb- With this starting 
point, one can apply the ideas of the proof of Theorem lA.il to show the existence 
of a solution of the equation (A -|- k'^)u = in R" such that ^{L) is a union of 
connected components of u~^{0), but it should be noticed that in this case the 
diffeomorphism <!> is not granted to be close to the identity. 

To illustrate Theorem lA.il let us consider the following easy application: 

Example A. 3. There is a function satisfying Au = u in R'^ whose zero set contains 
all the compact surfaces, that is, has a component of genus g for all nonnegative 
integers g (cf. Figure [5]) . This is an obvious consequence of Theorem lA.ll as there 
is a countable number of compact surfaces modulo diffeomorphism (indeed, they 
are customary labeled by their genus). 
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